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Abstract 

v 

During  this  period  I  have  extended  the  second-quantized  derivation  of 
quantum  i/f  noise  shown  in  the  Second  Annual  Report  to  the  general  case  of  N 
particles  present  in  the  final  state.  I  also  have  derived  the  quantum  i/f  cross 
correlations  and  the  corresponding  cross  -correlation  spectra,  which  are 
important  for  the  calculation  of  quantum  1/f  noise  in  kinetic  coefficients  such  as 
the  mobility  and  the  diffusion  coefficient  of  the  current  carriers  in  solids.  In 
order  to  better  explain  the  foundations  of  quantum  1/f  theory,  I  have  given  a 
derivation  of  the  quantum  1/f  Schroedinger  fields  from  quantum  electrodynamics 
with  the  use  of  coherent  states.  Finally,  I  have  given  a  direct  derivation  of  the 
quantum  1/f  effect  in  time  and  space.  In  terms  of  applications,  a  quantum  1/f 
noise  study  of  MIS  detectors  was  performed.  Experimentally,  with  the 
collaboration  of  the  group  of  Prof.  A.  van  der  Ziel,  an  excellent  experimental 
verification  of  quantum  1/f  theory  was  performed  on  semiconductor  diodes, 
transistors  and  vacuum  tubes,  and  a  review  article  on  the  results  of  the 
experimental  aplication  and  verification  of  my  theory  was  published  by  A.  van  der 
Ziel  in  the  Proceedings  of  IEEE  in  March  1938. 
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I.  INTRODUCTION 

The  present  report  will  give  a  general  description  of  the  very  extensive  work 
performed  in  the  reported  period.  This  work  was  done  along  two  major  directions: 

1.  Theory,  which  is  presented  in  Part  I  of  this  Report.  This  had  the  goal  to 
further  develop  the  quantum  1/f  theory  and  to  clarify  its  foundations  and 
starting  points.  This  part  has  nine  sections,  an  Appendix,  and  its  own 
references. 

2.  Application  to  MIS  infrared  detector  structures,  presented  in  Part  II.  This 
part  is  dedicated  to  a  detailed  quantum  1/f  noise  study  of  an  important  type  of 
infrared  detector,  the  metal  -  insulator  -  semiconductor  structure. 

Part  III  lists  the  publications  corresponding  to  this  grant  period. 

At  the  III  Conference  on  Quantum  1/f  Noise  and  1/f  Noise  in  Minneapolis, 


April  28-29,  1938,  Prof.  C.M.  Van  Vliet  also  presented  a  quantum  1/f  noise 
derivation  in  the  Van  Hove  limit.  At  present  our  attention  is  focussed  on 
applications  of  the  newly  calculated  cross  -  spectra  of  quantum  1/f  noise. 
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PART  I 


THEORY 


II.  DERIVATION  OF  THE  PAIR-CORRELATION  FUNCTION  FOR  BOSONS 
Denoting  by  <fr(x,t)  and  x(^.t)  the  single-particle  wave  functions  of  two 
scattered  bosons,  which  take  Into  account  both  the  interaction  with  the 
scattering  center  of  forces  and  the  simultaneous  interaction  with  the  photons, 
the  quantum  state  of  two  bosons,  both  emerging  from  the  same  interaction  process 
can  be  written  in  the  Heisenberg  picture  in  the  form 

|S>  *  (1/2)1/2  jd35jd3n^(r»0)x(n*0)'(»+(c)'|»+(n)  |0>,  (2.1) 

where  y(Z)  designates  the  Heisenberg  field  operator  at  a  fixed  time  t  *  0  and 
the  single-particle  wave  functions  have  also  been  taken  at  t  »  0.  The  two 
single-particle  wave  functions  Included  differ  only  through  some  external  and 
Internal  phase  factors.  Vectors  are  printed  In  boldface.  The  state  in  Eq. 

(2.1)  is  written  as  a  product  of  two  single-particle  wave  functions,  because  we 
neglect  all  Interactions  between  the  particles  except  for  the  quantum  exchange 
symmetry  between  Identical  particles.  We  have  limited  ourselves  here  to  two 
outgoing  particles  for  the  sake  of  simplicity  only;  the  general  N-partlcle  case 
Is  treated  In  Sec.  VIII. 

The  operator  of  the  equal-time  pair-correlation®  in  space  is 

0  «  tji+(xi)ij»+(x2)'l'(x2),J'{xl).  (2*2) 

This  corresponds  to  a  density  autocorrelation  function®.  The  presence  of  two- 
particle  coordinates  in  the  operator  0  does  not  mean  that  we  are  considering 
two-particle  Interactions,  It  only  means  that  the  expectation  value  which  we  are 
calculating  depends  on  the  relative  position  of  the  particles.  Using  the  well 


known  commutation  relations 

<Kx) <*»+(y)  -  ♦+(y)«i»(x)  -  6(x  -  y). 

(2.3) 

♦(x)ij»(y)  -  «|»(y)v(x)  -  0, 

(2.3a) 

♦+(x)/(y)  -  ♦+(y )t|»+(x)  -  0, 

(2.3b) 

we  obtain  the  matrix  element: 


<S°|0|S°> 

-  (l/2)<0|4.(n,)*(l,)*+(x1)/{x2)^(jr2)4-(x1)1i,+(|)4l+(n)|0> 

*  (1/2) [5(4'  -  xi)6(n‘  -  x2)  +  6(n'  -  3^ ) 6 ( C '  -  x2) ] 

[6(n  -  x2)<S(5  "  Xj)  +  6(5  -  x2)6(tT  -  Xj),  (2.4) 

where  |S°>  is  the  state  with  well  defined  particle  coordinates. 

The  pair-correlation  function  is  then 

A  =»  <S|0|S>  »  (l/2)<x*(x2)<>*(xi)^(^l)x(x2)  +  X^(x2)^*(xi)«(x2>x(^l) 

+  X*(x\)t*(x2)<>(xi)x(x2)  +  X*(xi)i»*(x2)l>(x2)x(xi)>.  (2.5) 

— _  > 

1)  If  we  assume  that  the  wave  functions  <j>(x),  x(x)>  ...of  different 
particles  in  the  outgoing  flux  differ  only  through  a  general  phase  factor,  we 
obtain  from  Eq.  (2.5) 

<A>  =*  2< 1 4>(xi )  |2  |<j>(x2)  |2>, 

which  is  similar  to  our  previous  result2  2<  J 4> (t ) <)> ( t+x )  |2>  a  2<  1  <t> ( t ^ ) <j> (t2 )  |2> , 
the  only  distinction  being  the  use  of  the  spatial  coordinate  along  the  beam 
instead  of  the  time  coordinate.  As  before*"^  we  write  from  Eq.  (2.5) 

<A>  =  2< 1 4>no ( x ! )  +  E-j  i>1Br(x1)  |2  Uno(*2)  +  $1(Jr(x2)  l2> 

-  2<( <t>no |2  +  li  U’Brl2)2 
+  ?1>*n0(xi)4>no(x2)  Ei  <4»*iBr(x2),l>1Br(xl)>  +  cc* 


(2.6) 


Here  we  have  separated  the  part  with  bremsstrahlung  into  the  photon  mode  1,  ^gr 
from  the  main  part,  <j>no,  without  bremsstrahlung,  and  cc  denotes  the  complex 
conjugate  of  the  preceedlng  term.  Each  ♦Sr  term  has  an  Independent  random 
phase;  as  we  shall  see  In  Sec.  IV,  this  Is  the  random  In  Itlal  phase  of  each 
photon  mode  1,  where  1  Includes  both  the  wave  vector  and  the  polarization.  The 
averaging  In  Eqs.  (2.6)  and  (2.6a)  Is  with  respect  to  these  phases,  or.  If  the 
emitted  photons  are  Included  Into  the  slngl e-partlcle  wave  function,  the  average 
Includes  also  the  expectation  value  In  the  space  of  photons.  Due  to  the  uniform 
motion  of  the  scattered  particles,  the  outgoing  radial  coordinate  and  the  time 
are  equivalent;  we  will  check  this  Intuitive  conjecture  later  In  Sec.  VI  and  we 
will  prove  by  direct  calculation  of  the  correlation  in  space  and  in  time  that  It 
Is  a  very  good  approximation. 

We  have  denoted  here  by  $  the  stochastic  Schroedinger  field  <j»  used  In  our 
previous  publications1"5.  This  field  was  not  second-quantized,  and  not  an 
operator  as  far  as  the  charged  particles  were  concerned.  It  was  always  called  a 
field  rather  than  a  wave  function.  In  order  to  emphasize  the  presence  of  random 
phase  factors  in  its  expression,  which  make  it  stochastic  in  nature,  while  the 
term  wave  function  was  reserved  for  pure  states  rather  than  mixture  of  states. 
The  stochastic  Schroedinger  field  description  is  equivalent  with  the  density 
matrix  description  of  quantum  mixtures.  The  reason  why  still  was  an  operator 
In  the  first  paper1  Is  that  while  the  particles  were  not  second-quantized  in 
that  paper,  the  electromagnetic  field  was.  Therefore,  in  that  paper  ij>  was  an 
operator  only  in  the  space  of  the  photons. 

2)  If,  on  the  other  hand,  the  single-particle  wave  functions  <j>(x),  x(x)>  ••• 
also  differ  through  mutually  independent  sets  of  random  phases  in  their  energy 
loss  (<^Br»  x^Br*  •••)  parts  as  we  conclude  in  Sec.  IV,  the  first  and  the  last 
terms  In  Eq.  (2.5)  are  constant  and  do  not  yield  1/f  noise,  but  the  two  middle 


f.l  *.f 
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terms  give  again  the  same  result  as  In  Eq.  (2.6) 

<A>  -  <U(xx)  |2>< ! ♦(x2)  |2>  +  !<<j>*(x1)4>(x2)>  |2 

“  (Unol2+  Z1  I <J»iBr I2) 2  +  l$*no(xl)<>no(x2)  +  Ef  V^Br^lJ^BrO^)  l>  I2 

*  (Unol2  +  Ei  I ♦1  Br I2) 2  +  Unol4  +  <21  I^Brl2)2 
+  /no(xl)<»no(x2)  Z1  Br( x2 ) t 1 B r< X1 )  +  cc 

*  < I ♦( xi )  |2><  | <t>(x2)  i2>  +  <U(xi)  |2  U(x2)  |2>.  (2.7) 

Both  In  the  first  and  last  form  we  notice  that  the  first  term  Is  constant,  while 
the  second  is  the  familiar  APSPSWF.  Therefore,  this  is  also  similar  to  the 
APSPSWF  result  obtained  In  our  previous  papers1"5,  but  this  time  the  fractional 
spectrum  will  be  two  times  smallerfdue  (see  Sec.  VIII)  to  a  1/N  factor  with  N*2; 
fermions  have  1/(N-1).  Here  the  average  Is  with  respect  to  the  sets  of  random 
phases  present  in  the  energy  loss  (bremsstrahlung)  contributions,  or,  if  we 
Include  the  photons  in  second  quantization  into  the  final  state  (see  Eq.  4.15), 
the  average  sign  in  Eq.  (2.7)  also  includes  the  expectation  value  in  the  Hilbert 
space  of  the  photons,  which  yields  exactly  the  same  result.  The  arguments  have 
been  omitted  in  some  of  the  terms  which  do  not  actually  depend  on  them  according 
to  Eq.  (4.15)  of  Sec.  IV..  Sec.  VIII  gives  the  APSPSWF  result  for  any  N. 

—9  -♦ 

We  conclude  that  both  if  the  outgoing  particle  wave  functions  $(x),  xU)» 
...»  differ  only  by  an  arbitrary  general  phase  factor,  and  if  they  also  have 
Independent  sets  of  phases  in  their  bremsstrahlung  energy  loss  parts,  our 
previous  method  of  calculating  the  autocorrelation  function  is  justified  on  the 
basis  of  the  present  second-quantized  calculation;  see  Sec.  IV-V  for  details. 

In  Sec.  IV,  we  will  determine  the  correct  single-particle  wave  functions  (or 
Schroedlnger  fields)  which  lead  to  Eq.  (2.7)  and  to  the  fractional  spectrum 
2aA/fN  for  N  outgoing  particles  (Sec.  VIII)  with  and  without  the  Inclusion  of 


-9- 


the  photons  Into  the  final  state  of  the  scattered  particles. 

Why  do  we  allow  for  so  many  possibilities  ?  Do  we  need  to  consider  both 
cases  with  and  without  the  photons  Included  Into  the  final  state?  Often  in 
electrophysics  one  describes  the  quantum  motion  of  an  electron  In  an  external 
time-dependent,  electromagnetic  field  of  force.  In  this  most  convenient 
description^  energy  Is  not  conserved,  and  the  electromagnetic  field  energy  Is 
not  Included  Into  the  Hamiltonian  which  is  time-dependent.  The  state  does  not 
include  the  field  either,  and  Is  a  nonstationary  electronic  state.  This 
corresponds  to  the  treatment^  presented  in  this  section,  i.e.,  without  the 
photons  Included  In  the  state  (Sec.  IV).  An  equivalent,  more  fundamental, 
approach  Includes  both  the  charged  particles  and  the  field  energies  with  their 
Interaction  Into  the  Hamiltonian.  Then  the  photons  must  be  present  In  the 
state,  as  we  did  before*.  The  two  methods  are  equivalent,  and  give  the  same 
results,  with  no  contradictions.  We  consider  here  both  methods. 

Finally,  why  do  we  need  to  consider  both  the  case  of  identical  single- 

particle  wave  functions  (except  for  a  general  phase  factor)  and  the  case  of  wave 

functions  which  have  different  sets  of  phases  in  their  energy-loss  parts? 

Because  both  cases  may  occur  in  practice.  Indeed,  usually  the  incoming 

particles  In  a  scattering  experiment  have  a  random  shift  in  space  or  time  which 

scrambles  their  energy-loss  phases  and  eliminates  cross  terms  as  if  these  sets 

of  random  phases  would  be  completely  independent  for  different  particles. 

However,  in  emission  processes,  such  as  o-decay,  the  single-particle  wave 

functions  are  not  shifted;  they  all  start  leaking  out  at  t  *  0  when  the 

radioactive  source  was  prepared  by  chemical  separation.  Therefore  we  must 

consider  all  cases  in  spite  of  the  Inconvenience.  For  instance,  by  writing 

explicitly  <p(x)  *  9no(x)  +  z  ?iBr(x)  and  XU)  ■  Xno(x)  +  £  xiBr(x)a+i»  we 

1  i 

obtain  from  Eq.  (2.5)  in  which  the  asterisks  are  interpreted  as  crosses 
Indicating  hermitlan  conjugation,  the  same  result  as  in  Eq.  (2.7),  by 


.  I  ,  !i  ■! 


tv’* 
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Interpreting  the  brackets  as  vacuum  expectation  values. 


III.  DERIVATION  OF  THE  PAIR-CORRELATION  FUNCTION  FOR  FERMIONS 


In  the  case  of  fermions,  the  state  of  two  scattered  particles  emerging  from 


the  same  Interaction  Is  similar  to  Eq.  (2.1) 

| Sss •  >  -  (1/2)1/2  Jd3Cjd3^(C,0)x(n.0)/s'(t)>|»+s'(n)|0>,  (3.1) 

but  now  the  field  operators  satisfy  anticommutation  relations: 

*s(*)tV(y)  +  ♦V(?)*s(*)  s  6(*  ■  y)«s,s'»  (3*2) 

♦s(x)*s'(y)  +  *s*  (?)*$(*)  “  °»  (3.2a! 

♦+s(x)*V(y)  +  t+s*(y)*+s(x)  3  o.  (3.2b; 


(3.2a) 


(3.2b) 


The  operator  of  the  pair-correlation  function  Is  now  written  In  the  form 
0  -  Ess*  *+s  (x‘i)l»+s,(x2)’l's,(x2Hs(xl)*  (3.3) 

Its  expectation  value  In  the  state  given  by  Eq.  (3.1)  is  calculated  with  the 


help  of  the  matrix  elements 


<stmo> 

3  <  0 1  tt(  n ' ' )  *t(  ?  )  *t+  ( x  1 )  tt+  ( x  2 )  %( x  2 )  *,( x  1 )  V* (T )  %+  ( n )  1 0  > 
3  C-«(n'  -lq)5(E'  -  x*2)  +  6^'  “  xl)5(n‘  -  xg)] 
C«Cn  -  x*2 )«(V  -  xi)  -  6(T  -  *X2)«(n  -  xi)3; 

3  <  0 1  <j>  ( n ' )  <t»  ( C ' )  n>.+ ( x  l )  ( x  2 )  ( x  2 )  ♦.( x  1 )  U )  ^  ( n )  1 0 > 

|  f  _ ^  r  ^  v  r  T  w 

-  6("n'  -  ~X2) 5( C*  -  xl)«U  -  X2)«U  -  xl); 

<S°J0|S!> 

*r  n  it 

3  <0|*(n,)*(e,)t?(*l) ¥+( x2) ♦,( x2 ) ♦.( X1 ) ) <iv+( n)  |0> 

1  '  *  '  *  T  *  r 

3  fi(c‘  -  X2) 6 (n 1  -  xi)sU  "  x2)5(n  ”  xl)* 


(3.4) 


(3.5) 


(3.6) 


We  also  obtain  three  similar  expectation  values  with  all  spins  reversed. 


The  spin-averaged  pair-correlation  function  is  then 


A  3  (l/4)ls#S*  <Sss'|0f+  y  0|J+  <|JSss.> 


cv;*- 


*  (l/4)[|x(x\)  I2I^»(x2)  I2  ■  X*(xl)$*(x2)x(x2)<l,(xl) 

-X*(x2)^*(xl)x(xl)|<'(x2)  +  lx(xl)  I2  i*(x2)  I2 

♦  Ix(xJ) l2l*(*xi) I2  ♦  IxUi)I2U(*2)I23 

-  (1/2)Cx*(x2>^^(xi )«t»(xi)x(x2)  +  x*(xl)$*(x2Mx2)x(xl>J 
-  ( 1/4) Cx*(xi) <j»*(x2)4>(xi) x(x2)  +  x*(x2H*(xlMx2)x(xl)]»  (3.7) 

1)  If  we  assume  that  the  wave  functions  $(x)  and  x(x)  differ  only  through  a 
general  phase  factor,  we  obtain  from  Eq.  (3.7) 

<A>  -  (1/2)<|*(71)|2|*6£)|2>  (3.7a) 

which  corresponds  to  our  earlier  APSPSWF  result2  in  the  time  domain,  and  which 
yields  again  the  fractional  spectrum  2aA/f  found  in  our  previous  papers.  (See 
Sec.  IV  for  the  definition  of  aA.) 

2)  If,  on  the  other  hand,  the  single-particle  wave  functions  also  contain 
mutually  Independent  sets  of  random  phases  In  their  energy  loss  parts,  we  obtain 
from  Eq.  (3.7) 


<A>  =»  <|$|2><|$|2>  -  (1/2)  |<^*(xi)0(x2)>l2 


3  (Unol2  +  Z1  I^Brl2)2  “  (1/2)1  Uno I4  +  (Ei  UVl2)2 
+  t*no(*l)Wx2)  si  /iBr(x2)^iBr(xl)  +  cc) 

*  <|4>|2>2  -  (1/2)<|*(m)I2|*(^)I2>.  (3.7b) 


Both  in  the  first  and  the  last  form  above  we  notice  that  the  firstterm  is 
constant,  while  the  second  is  the  familiar  APSPSWF.  Therefore,  this  result  is 
again  similar  to  our  previous  APSPSWF  result*-2  shown  in  Eq.  (2.6)  or  (3.7a), 


*»•*•»*'  u* 
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and  yields  the  same  fractional  spectrum  2aA/f.  However,  at  small  distances  we 
now  obtain  something  similar  to  a  Fermi  hole  due  to  the  minus  sign  in  Eq. 

(3.7b).  The  minus  sign  corresponds  to  a  180°  phase  shift  or  delay.  This  means 
that  the  fermions  avoid  each  other  at  close  range.  Over  large  distances,  l.e., 
at  the  low  wave  numbers  significant  for  1/f  noise,  the  distribution  and  the 
spectral  density  are  not  affected.  If  we  explictly  Include  the  photons  into  the 
final  state  and  multiply  ^Br  by  a  photon  creation  operation  a^j,  we  obtain 
again  exactly  the  same  result  shown  In  Eq.  (3.7b).  This  shows  that  although  the 
phase  shift  caused  by  the  Pauli  principle  at  short  distances  was  not  Included, 
our  previous  treatment*"5  was  correct  because  It  started  with  exactly  the  same 
product  of  single-particle  wave  functions  and  therefore  yielded  the  same  final 
result  as  the  present  calculation.  Whether  this  result  should  simply  be 
interpreted  In  terms  of  frequency  beats  Is  a  philosophical  question;  I  know  It 
should.  However,  while  the  beats  are  caused  by  exactly  the  same  Interference  of 
each  single-particle  wave  function  with  itself  which  we  had  considered  from  the 
beginning,  they  occur  because  the  single-particle  wave  functions  are  Included  as 
symmetrized  products  into  the  many-particle  wave  function  of  the  outgoing  flux. 
The  Fermi  hole  is  absent  In  Eq.  (3),  because  there  only  the  antiparallel  spins 
contribute. 

We  finally  conclude  for  fermions  as  well  as  for  bosons,  that  both  if  the 
outgoing  particle  wave  functions  <j>(x),  *(x),  ...»  differ  only  by  an  arbitrary 
general  phase  factor  and  If  they  also  have  independent  sets  of  phases  in  their 
bremsstrahlung  energy  loss  parts,  our  previous  method  of  calculating  the 
autocorrelation  function  is  justified  on  the  basis  of  the  present  second- 
quantized  calculation;  see  Sec.  IV-V  for  details.  All  claims  to  the  contrary 
which  have  been  recently  voiced,  are  therefore  wrong.  We  now  need  to  insert  the 
correct  single-particle  wave  functions  explicitly  into  Eqs.  (2.5)  and  (3.7). 
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Before  we  do  this  direct  calculation,  we  briefly  discuss  the  connection 
between  density  and  current  density  autocorrelations.  The  current  density 


operator  in  second  quantization  is 


J  «  Crt/21mH+(J)($-$)*(x), 


(5.6) 


where  now  ♦(x)  and  *|»+(x)  are  particle  field  operators,  while  the  same  notation 
was  used  earlier1  to  designate  single  particle  wave  functions  which  are 
operators  in  the  Hilbert  space  of  photons,  because  the  current-carrying 
particles  were  not  second-quantized,  while  the  emitted  photons  present  In  the 
final  state  were  quantized.  The  operator  of  the  current  correlation  which 
replaces  Eq.  (2.2)  Is 


Oj  *  (V2m)2Y+(xiH+(x2)($2‘^2)|i’(*2)(^l“^l)'l'(*l)» 
while  In  Eqs.  (2.5)  and  (3.7)  we  have  to  make  the  substitutions 

♦  (*!)'-  fl/21«)(MlM*l>  •  (**/n»M*l) 


(5.7) 


♦(*2)  ♦  (V21m)($2-$2H(*2)  *  (15K/m)*(St2)  (5.8) 

and  similar  for  XU),  while  and  x*  are  left  formally  unaffected.  Here  we 
have  used  the  form  of  $(x)  and  x(*)  given  by  Eq.  (5.1)  or  (5.2)  which  allow  only 
for  very  small  momentum  deviations  q  «  K  caused  by  infraquanta,  which  can  be 
always  neglected.  Therefore,  Eqs.  (2.5)  and  (3.7)  will  only  acquire  a  general 
factor  (TiK/m ) 2 .  We  conclude  that  in  general 

<S|0j|S>  -  (TiK/m)2<S|0|S>,  (5.9) 

l.e.,  that  the  current  density  autocorrelation  Is  proportional  to  the  density 
autocorrelation  If  the  single-particle  states  are  close  to  momentum  eigenstates. 
The  corresponding  spectral  densities  are  also  proportional,  with  the  same 
proportionality  factor  (HK/m)2.  The  fractional  autocorrelations,  however  are 
the  same  for  the  current,  as  for  the  density,  because  they  have  the  squared 
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average  value  In  the  denominator.  The  same  equality  holds  for  the  corresponding 
spectra. 

Finally,  for  single-particle  states  which  are  superpositions  of  states  close 
to  momentum  eigenstates,  the  above  proportionality  remains  approximately  valid, 
if  we  replace  the  proportionality  factor  with  <(fiK/m)2>,  where  the  average  Is 
over  the  distribution  of  K  . 
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IV.  DERIVATION  OF  THE  SINGLE-PARTICLE  WAVE  FUNCTIONS 

The  derivation  of  the  Schroedinger  field  will  be  performed  with  the 
help  of  the  Green's  function  method  similar  to  the  method  used  in  a  earlier 
calculation  by  Kroll  and  Watson7,  extended  to  the  case  of  interaction  with 
all  electromagnetic  modes  of  the  universe*. 

It  is  most  convenient  to  describe  the  electromagnetic  field  in  terms 
of  plane  waves.  The  vector  potential  is  taken  in  the  radiation  gauge  as 

A(r,t)  r  2  ♦  a*k,i(t)e-“-r].  (4.1) 

The  polarization  vectors  ut,i  and  Uk,a  are  mutually  orthogonal  unit  vectors 
perpendicular  to  k. 

Hie  Schroedinger  equation  for  an  electron  moving  in  a  vector  poten¬ 
tial  A  and  scattering  potential  V  is 

( l/2m)  [-ifiv  -  eA/c]  it/>  +  Vtf  =  ifitf.  (4.2) 

A  dot  has  been  used  to  indicate  the  time  derivative.  The  electromagnetic 
field  is  treated  as  a  classical  field  at  this  point.  In  order  to  eliminate 
the  A*  term  from  Eq.  (4.2),  we  write 

V  -  expn-i/hWMeV^OA’dt’]#  . 

Thus,  Eq.  (4.2)  is  reduced  to 

[(-fiV2m)V*  +  ( ieft/mc)A*V  +  V]«t>  =  ih<t>  .  (4.3) 


It  is  convenient  to  consider  first  the  influence  of  a  single  electromagnetic 
mode,  i.e.  a  single  term  from  Eq.  (4.1).  Therefore,  we  take  A  »  acos ( ctft+7) , 
where  y  is  an  initial  phase  constant,  and  we  treat  V<X>  as  a  perturbation 


source  tens.  The  solution  for  Eq.  (4.2)  is  an  incoming  plane  wave  plus 
scattered  waves,  given  by  the  integral  equation 


<J>kJr,t)  -  4>kc  -  XdJx’.p  dt'GVfr’)#  kc(r\t*).  (4.4) 

Here  is  the  solution  of  the  homogeneous  equation,  i.e.  with  V  =  0,  and 
can  be  written  in  the  form 

tfke  s  e»»efexp[-(ifi/2m)J*(k*,  -  2ek0- A/ch)dt] .  (4.5) 

G  is  the  Green's  function  which  satisfies  the  equation 

[  (-ha/2m)V  *  +  (iehA-V/mc)  -  ih<5/0  t]G  =  $(r-r*  ) .  (4.6) 


Given  A  =  acos(<ut  +  y) ,  G  can  be  found  to  be 

I 


| 


G  s  [i/(2x)  ^JJd^e1  *(r-r'  )exp(-(ifi/2m)  [kat  -  2ek*asin(<ut  +  7)/hca>]} 
x  exp[(ih/2m)(kat'  -  2ek-asin(crt(  +  r)/Rcft>],  (4.7) 

In  the  first  Bom  approximation  we  set  <D  ^  ( r ’ , t ’ )  =  4>\o  (r’  ,t’ )  in  the  in¬ 
tegral  present  in  Eq.  (4.4)  and  obtain  for  the  scattered  wave 


=  [i/(2Ji)*ft]Jdax’J»  dt,V(r’  )Jdake*  r“r ’  >exp{~(iii/2m) 

[kat  -  2ek*asin(o>t  +  7)/hc©]}  x  exp[  (ifi/2m)  (ha  -  kao)t’] 
x  {exp(ie(k«  -  k)asin(ot'  +  r)/2lncal]e,  kc'r ’)  • 

Using  the  relation 

=  £  J»(/9)e“»<*‘’*r>f  (4.8) 

where  J«(/3)  is  the  nth  order  Bessel  function,  we  expand  the  expression  con¬ 
tained  in  curly  brackets  in  Fourier  series.  Then  Eq.  (4.4)  takes  the  form 

(r,t)  -  dk8(r,t)  =  [-i/(2x)aft]Jdax,_J‘dt'V(r')Jdake‘fc(r-r*, 

x  exp{-(ifi/2ra)[kat  -  2ekasin(«t+r)/hca)J) 

x  exp((ifi/2m)(ka  -  ka)t’][£  J«09)e*  ■f*'1 '  ♦  ^]e‘  ’ .  (4.9) 

After  performing  the  integration  over  t’  we  use  a  contour  integration  method 
for  k.  Then  Eq.  (4.9)  is  reduced  to 

<r,t)  -  < r,t)  =  (-m/(2x)ha]Z  JB(/3)e‘«*[e‘*<»  >'/r] 

e  °  n«-*» 

x  exp(ih[ka(n)t-2ek(n)asin(cjt  +r)/hc<a]/2m)  x  Jdax'e 1  » > r  ’V(r*  )el  k«  r  * , 


where 


$  -  -e(tik*  -  tik o )  •  a/mcfia)  =  -e<J-a/mchG>,  (4.11) 

and  Q  is  the  momentum  transfer.  In  Eq.  (4.10)  k(n)  is  defined  by 

(Tvk(n) )  */2m  =  (hko)*/2m  -  k.  =  k(n)  =  k(n)r/r.  (4.12) 

The  total  scattered  wave  can  be  written  as 

0*  =  [-m/(2jr)f»Jr]2  rexp{-ililk2(n)t  -2ek(n)  •asin(wt+r)/ftc<u]/2m} 

x  Vk(«)  ,k  J.(y3)e‘«r,  (4.13) 

where  Vk(«),k  =  J*e* 1  k(  ■> r  *V(r'  )e*  *o*  r  *d1x’  is  the  scattering  matrix  ele¬ 

ment  calculated  without  consideration  of  the  interaction  with  the 
electromagnetic  field  oscillators.  Gr*ertliittic*.  to dl  fold aullldc riSfiVe<  */>**«  /.&re*cJv, 
So  far  the  electromagnetic  field  has  not  been  quantized  and  was  con¬ 
sidered  as  a  classical  field.  We  are  interested  in  the  corresponding  ex¬ 
pression  of  the  scattered  single-particle  wave  function  when  the  electromag¬ 
netic  oscillators  are  quantized.  Therefore  we  first  linearize  Eq.  (4.13) 
with  respect  to  the  electromagnetic  potential  wherever  a  dependence  on  g  is 
present: 

0*  =  [-m/(2x)ttIrJVi,k  ei*r-,s‘''*exptieK'asin(<ut  +  rJ/cmu] 

{1  +  el  2} 


=  [-m/(2x)fi,r]Vf ,k  elIr-lE*/*exp[l  +  eiK-asin(<at  +  r)/cm6>] 

{1  +  e1  (iyt**r+<‘>eQ*a/2mCficu  -  e*‘ <<ut*^r4^‘)eQ-a/2mch6)} .  (4.14) 

Since  k  =  (v/c)q<<q,  the  r-depedence  of  the  el'ctromagnetic  potential  can  be 
neglected,  as  it  is,  in  exp  [ieK*asin(<ut  +  r)/cm<u]. 


Here  we  have  introduced  the  notations  E  =  (fiK)*/2m  and  K  =  k(0).  In  the 

last  form  we  will  neglect  the  terra  with  sin(<at  +  r)  tfcat  corresponds 
part  cftA* 

to  a\coherent  quantum  1/f  noise  contribution  which  has  been  considered  else¬ 


where  before"1 12  , 


because  this  term  can  be  considered  constant  of  negli¬ 


gible  magnitude  in  the  calculation  of  equal  time  spatial  correlations.  We 
also  neglected  the  small  difference  between  k(l)  =  K  -  afi/vK  and  K  in  $  and 
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vk ( 1 ) ,ka •  We  conclude  that  the  quantization  of  the  electromagnetic  field 
transforms  into  an  operator  in  the  photon  Hilbert  space 

^(r,t)  =  (C/r)e,“r-,,‘''#{l  -  2  k,  ib*(k,l)expi(-4rt:  +  qr  -  7i)ak,i 

+  2  k,  ib(k,l)expi(<ut  -  qr  +  7i)a*k,i).  (4.15) 

This  is  the  form  of  the  single-particle  wave  functions  which  are  used  in  the 
present  paper  with  or  without  the  photons  included.  Here  we  have  introduced 
b(k,l)  ■  (1/2)8  and  the  constant  C  which  designates  the  factor  In  front,  and  a 
sum  which  Includes  all  electromagnetic  modes  with  annihilation  operators  a£ti« 
In  the  space  of  the  electron  states  $(r,t)  is  just  a  single-particle  wave 
function.  We  have  denoted  by  q  the  small  decrease  In  the  particle  momentum 
required  by  Eq.  (4.12).  We  have  q  3  (K/E)*Hu  3  ck/v  *  w/v,  with 

(L/2*)J-4*-2  i<|  b(k,l)  |  *>k*dk  =  e*Q*aV4m*c*!\*<iJ1-2kIdca/3c  =  <xAda>/a,  (4.16) 


where  a  =  e*/4rtfic  =  1/137,  A  =  (2Ql/3«nJc*) ,  and  where  <  >  is  an  angular 
average.  We  have  considered  the  spontaneous  emission  caused  by  vacuum  fluc¬ 
tuations  only,  yielding  far  the  amplitude 


(4.17) 


The  annihilation  part  is  included  in  Eq.  (4.15),  but  does  not  contribute  to 
the  quantum  1/f  noise  on  the  background  of  the  electromagnetic  vacuum.  If 
the  calculation  is  performed  on  the  thermal  radiation  background,  however, 
we  get  a  white  noise  contribution  added  to  the  quantum  1/f  noise  which 
remains  the  same*.  Here  we  have  performed  the  transition  from  just  one 
electromagnetic  mode  to  the  general  case  with  all  electromagnetic  modes  ad- 
hoc,  but  in  our  previous  paper*  this  transition  was  presented  in  detail. 


In  Eq.  (4.15)  we  notice  the  presence  of  the  random  phases  n  which 
where  introduced  as  initial  phases  of  the  electromagnetic  oscillators,  and 


which  are  independent  for  each  electromagnetic  mode  i  (of  given  polarization 
1  and  wave  vector  k)  of  the  universe.  Since  the  various  scattered  particles 
are  independent  of  each  other  and  of  the  electromagnetic  modes,  we  have  to 
consider  the  set  of  random  phases  T\  different  and  independent  in  the  wave 
function  of  each  particle.  tke  Jiscu-ssCcn.  befcce  ifttr  fg, (A/.u), 

V.  CALCULATION  OF  THE  PAIR-CORRELATION  FUNCTIONS  AND 
OF  THE  CURRENT  AUTOCORRELATION 


In  Sec.  IV  we  have  determined  the  single-particle  wave  functions 
in  the  form,  valid  at  t  =  0, 

<t>(x, t)  =  (C/x)exp(iKx)  [1  +  2  k ,  ib(k,l)exp(-iqx)  ] ,  (5.1) 

or,  if  the  photons  emitted  in  the  bremsstrahlung  process  are  included  in 
second  quantization, 

tf(x,t)  =  (C/x)exp(iKx) [1  +  2  a, ib(kj)exp(-iqx)a*k, i],  (5.2) 

Substituting  these  expressions  into  the  calculated  expectation  values  in  Eq. 
(3.7),  we  obtain  for  the  case  of  fermions 

A(x»,  x2)  =  (C*C/xiXj) *{1/2  +  2  a ,  i  |  b(k,l)  |  *[2-cosq(xx  -Xj)] 

+  2  kkMi  •  |  b(k,l)  |  *  |  blk’,1’)  |  *[1  -( l/2)cos(q  -q’ )  (xi  -  x2 )  ] } .  (5.3) 

Note  that  if  Eq.  (5.2)  is  substituted,  the  expectation  on  the  electromag- 
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net  ic  vacuum  state  is  calculated.  If  the  thermal  radiation  background  is 
also  taken  into  account,  a  small  white  noise  term  is  obtained  in  addition  to 
the  1/f  noise*.  The  constant  part  is  the  squared  expectation  value  of  the 
concentration  of  particles.  Dividing  the  variable  part  by  this  constant 
part, we  obtain  the  fractional  spectral  density  of  the  particle  concentration 
n,  or  current  density  j  (with  LJ  being  the  volume  of  the  normalization  box 
considered) 

S.(k)dk/n*  =  Sj(k)dk/j*  =  22  i<<|  b(k,l)  |  *>>(L/2*)3  4xk*dk 
/{I  +  42  k.i  |  b(k,l)  |  *  +  22  kk*n»|  b(k,l)  |  *|  b(k*,l»)  I  *) 

=  22  i«|  b(k,l)J  *>> (k*/2wJ )L3dk  =  2aAdk/k  =  Sj(f)df/j*  =  2aAdf/f,  (5.4) 

which  is  in  agreement  with  our  previous  results,  and  which  also  includes  a 
180*  phase  shift  due  to  the  exclusion  principle  which  is  important  only  at 
short  distances  between  the  particles.  In  the  final  form  we  have  trans¬ 
formed  to  the  frequency  f. 

The  expression  of  b(k,l)  used  in  Eq.  (5.4)  was  derived  in  Sec.  IV. 

It  is  the  well  known  expression  of  the  bremsstrahlung  amplitude  in  any  scat¬ 
tering  process 

b(k,l)  =  -  (e/n«iu)(h/a.3)w*ApAi,  (5.4a) 

where  Ai  is  the  polarization  vector  of  the  mode,  and  Ap  is  the  momentum 
change  of  the  particles  in  the  scattering  process  considered.  In  Eq.  (5.4) 
an  average  <<  >>  over  the  angular  variables  in  k  and  a  summation  over 

photon  polarizations  were  performed,  leading  to  the  usual  definition  of  the 
quantum  1/f  noise  coefficient,  also  known  as  the  nonrelativistic  form  of  the 
infrared  exponent  defined  in  quantum  electrodynamics 

aA  -  (2o/3x) (A  p/mc) *,  (5.4b) 

with  a  being  Somnerf eld's  fine  structure  constant  a  =  e*/ftc  =  1/137. 


In  the  case  of  bosons  we  substitute  Eq.  (5.1)  or  (5.2)  into  Eq. 

(2.5b)  and  obtain 

A(xi.xj)  =  2  +  22  k.  i  )  b(k,l)  |  *£!-»■  cosq(xi  -  x»)] 

♦  2  *,i2  »•.!•!  b(k,l)  |  *|  b(k’,l')  |  *[1  +  cos(q  -  q’)(xi  -  x*)]; 

Sa(k)dk/(n) 1  =  Sj(k)dk/(  j) 2  =  2  i«|  b(k,l)  |  *>>(L/2 *) »  4sk*dk 

=  aAdk/k  =  aAdf/f.  (5.5) 

Here  at  short  distances  we  notice  an  increase  of  A(xi,  xi).  The  fractional 
spectral  density  is  reduced  by  factor  of  two  compared  to  the  case  of  fer¬ 
mions.  Both  results  generalize  our  previous  results  to  the  case  of  short 
distances  and  prove  the  existence  of  the  conventional  quantum  1/f  effect  in 
second  quantization,  even  if  the  photons  are  included  into  the  final  state 
(which  is  to  be  avoided  if  the  electromagnetic  field  is  included  in  the 
Hamiltonian  as  an  external  field,  as  we  did  in  earlier  papers) . 

The  transition  from  the  wave  number  spectrum  to  a  frequency  spectrum 
in  the  last  form  Eqs.  (5.4)  and  (5.5)  is  based  on  the  equivalence  of  the 
density  distribution  along  the  outgoing  rays  with  the  time-dependence  in  a 
fixed  point  along  the  flux  of  particles  which  all  move  with  the  speed  v. 

The  coordinate  x  along  the  ray  and  the  time  t  are  in  this  case  equivalent 
variables,  and  therefore  the  wave  number  k  and  the  frequency  f  are  also 
equivalent  variables  for  the  spectral  density.  We  can  therefore  write  2xf  = 
vq  =  ck  and  S(k)dk  =  S(f)df,  which  justifies  the  last  form  of  Eqs.  (5.4)  and 
(5.5).  In  order  to  check  the  validity  of  this  procedure,  we  will  now  per¬ 
form  a  direct  calculation  of  the  pair-correlation  function  in  time  and 


space. 


VI.  DERIVATION  OF  THE  PAIR-CORRELATION  FUNCTION  IN  TIME  AND  SPACE 


Working  in  the  Heisenberg  representation  as  above,  we  can  generalize 
the  derivation  presented  above  by  including  two  different  tiroes  in  the 
operator  of  the  pair-correlation,  although  this  requires  more  calculation. 

To  simplify  the  integrals,  we  consider  again  the  state  of  two  outgoing  fer¬ 
mions  very  far  from  the  place  where  they  have  independently  suffered  the 
same  interaction,  so  that  the  outgoing  spherical  waves  can  be  approximated 
by  plane  waves.  Starting  therefore  with  plane  waves  similar  to  the  spheri¬ 
cal  waves  used  in  Eq.  (3.1),  with  q  =  qK/K  =  ckK/Kv  and  q’  =  ck’K/Kv,  we  ob¬ 
tain  again  the  pair-correlation  function  in  the  form 

A(xi ,  ti;  xj,  t,)  =  l/2Id»$d3*iJ£’dV 
{exp(-i&p’)  +  2  k  • , »  ‘ySMk’  ,s’  )exp[-i(K-q* )rf  ]a* • ,  •) 

(exp(-iK£*)  +  2  k,ib*(k,s)exp[-i(K-q)5']ak.} 

(exp(iKg)  +  2  k.  •b(k,s)exp[i(K-q)^)a*k*} 

{exp(iK^)  +  2  k*,f^(k,,s,)exp[i(K-q')7]a*k**>) 

(1/4)2  S3 *<S°ss •  |  Of|+  <y  041+  04|  ,t  i  ,2.i2  |  S*ss *  >.  (6.1) 

The  operator  of  the  pair-correlation  function  now  contains  two  consecutive 
times.  The  creation  and  annihilation  operators  for  particles  obey  anticom¬ 
mutation  relations  similar  to  Eqs.  (3.2)-(3.2b) .  Using  these  operators,  the 
field  operators  can  be  expanded  in  terms  of  plane  waves 

tfs(r,t)  =  V-*'*2  p  expi(pr-<wt)c»s,  0*s(r,t)  =  V-»/J2  »  expi(-pr+<ut)c**s(6.2) 


which  also  contain  their  time  dependence,  as  needed  for  the  operator  of  the 
pair  correlation  function  in  the  Heisenberg  representation.  This  operator 
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has  an  expectation  value  which  can  be  writtenf  for  spin  up  only,  in  the  form 
<S»M  |  Ofxt,  ti,  Xa,  ti)  |  S«ff> 

=  <0  |  0(-7'  )$(£’  )0*(xi,ti)0*(xa,ta)^(xi,tj)0(xi,ti )0*(S)0*(7)  I  0> 

=  (1/V«)2  ■■■•■•£  .,  expi[a*7'  +n'£'  -/<’xi  +<u(ju,)ti  -j/'xa  +o>(i/’)ta 
+/Kj  -<u(^i)ta  +1*1  -<u(i/)ti  -m7  -n£] 

<0  |  Cn>Cii'Cy*,  C^,  Cv  C*»C* a  |  0  > 

=  (M/h)«t1-Ha-Jexp{(iM/2fi)[(x1-f)Vti  +(xa-7)*/ta 
-(xi-7')*/ti  -(xa-^)Vta]} 

+(M/h)«ti-»ti-iexp{(iM/2fi)[(xi-7)*/ti  +(xa-|)*/ta 
-(xi-^')Vti  -(xa-7')Vta]} 

-(M/h)«ti-»ta-»exp{(iM/2K)[(xi-$) */ti  +(xa-7) */ta 
-(xi-e’)Vti  -(x,-7»)  Vta]} 


-(M/h)«ti-Ha'*exp{(iM/2K)[{xi-7)Vti  +(xa-f)*/ta 
-(xi-7’)Vti  -<x,-r)Vt,]>. 


(6.3) 


Substituting  into  Eq.  (6.1),  and  integrating  with  respect  to  7,  f  and 
7’ ,  we  obtain  for  the  part  with  spin  up  only 


8Ajj(xi,  ti,*xa,  ta) 

s  <0(1  +2  n«/9*(k’  )expi(q’xi  -1iq’(K  -q’/2)ti/M]a*») 

(1  +2  *b*(k)expi[qxa  -Dq(K  -q/2)ta/M]ak) 

(1  +2  »b(k)expi[qxi  +hq(K  -q/2)tj/M]a* k) 

(1  +2  k./3(k*)expi[q’xa  +tiq’-(K  -q’/2)ta/M]a»k«}0> 

=  1  i'|  b(k’)|  Jexpi(qt(Xi  -xa)  -ttq'-(K  -q’/2)(ti  -ta)/M] 

+  2  k  |  b(k)  |  Jexpi[q(xa  -Xi)  -tiq(K  -q/2)(ta  -ti)/M] 

+  2  k.k- |  b(k)b(k»)  |  »expi[(q  -q’Hxa  -Xi)  -fiK(q  -q’)(ta  -t»)/M 
+ti(q*  -q’»)(ta  -t,)/2M]. 


(6.4) 
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The  first  form  of  this  result  proves  the  presence,  with  ttt*v  of 
four  single-particle  functions  (APSP5WF  combination)  in  A^. 

In  a  similar  way  the  expectation  values  corresponding  to  Eqs.  (3.5)  and 

(3.6)  are  calculated  and  integrated  with  respect  to  £,  n,  €'  and  /?’ .  Put¬ 
ting  together  the  terms  for  all  spin  orientations  as  before,  we  finally  ob¬ 
tain  the  desired  pair  correlation  function  in  the  form 


A(xi,t»;x»,tj)  =  1/2  +  2  a , i  |  b(k,l)  |  *{2  -  cosq[xi-xa  -v'(ta-ta)]} 
+2  am ■  1 1  •  |  b(k,l)  |  *  |  b(k*  ,1)  |  Ml  -<l/2)cos(q-q*)-txi-xa  -v*  ’ (ti-ta) ]) ;  (6.5) 
v’=  v( 1  -q/ZK)  =  v( l-e/4E)  =  v;  v’ ’=  v[l  -(q+q*)/2K]  =  v[  1  ~(e+e’ )/4El  =  v. 


The  approximations  v*  =  v  and  v’  ’  =  v  are  justified  because  the  soft  photon 
energy  e  -  4- 10-1  seV  for  1Hz  is  negligible  compared  to  the  energy  E  of  the 
particles  which  may  be  of  the  order  of  leV.  This  result  directly  proves  our 
heuristic  generalization  of  Eq.  (5.4)  to  the  case  of  different  times  in  the 
pair-correlation  function.  The  extension  to  the  case  of  bosons  is  trivial. 


VII.  CK3SS -CORRELATIONS  AND  SPECTRA 


In  practical  calculations  of  quantum  1/f  noise  in  condensed  matter 
and  electronic  devices  we  often  need  the  cross-correlation  of  the  differen¬ 
tial  scattering  cross  sections  at  different  angles.  This  is  in  fact  the 
cross-correlation  of  the  outgoing  current  densities  scattered  into  different 
directions,  with  different  incoming  wave  vectors  Ki  and  Ka  considered.  If 
we  restrict  ourselves  again  to  spatial  correlations  for  simplicity,  we  need 
to  examine  the  cross-correlation  of  the  probability  densities  in  the  scat¬ 
tered  wave  at  various  distances  from  the  scattering  center  and  in  various 
directions . 


For  bosons  we  start  from  Eq.  (2.5)  and  consider  ^(x)  and  *(x)  as  in 


Eq.  (5.2),  but  different  from  eficb  ether; 

<£(x,t)  =  (C/x)exp(iKix)[l  +  Z  n, ib(k)exp(-iqix)a*k, i ] ,  (7.1) 

*(x,t)  =  (C/x)exp(iKjx)  [1  +  Z  k.  i/9(k)exp(-iq»x)a*k,  j  ] ,  (7.2) 

where  now  b  corresponds  to  the  momentun  change  K'  -  Ki,  and  $  corresponds  to 
the  momentum  change  K"  -  Kj  of  the  second  particle,  in  addition  of  having  a 
independent  random  phase  different  from  that  in  b,  for  each  k,  due  to  the 
different  initial  phases  of  the  electromagnetic  oscillators  registered  by 
the  second  particle.  As  in  Eq.  (2.5),  we  obtain  now 

A(Ki  ,Ki,K'  ,K"  ;xi  ,xa)  =  <S|0|S>=  ( 1/2)  {  |  *(xi)  |  *  |  tf(x,)  |  * 

♦  x*(xt)<t*(xi)x(xi)<t>(xi)+xHxt)<t*(xt)z(xi)<f>{x»)  +|  2(x«)  I  *|  4>(xi)  |  *} 

=  [1  +Z  k.i  I  /9(k,l)  |  *111  +  Z  k.i  |  b(k,l)  |  *] 

♦  (l/2)exp[iK* (xa-xi)  +  iK"(xi-Xi)][l  +  Zk.il  ,fl(k,l)  |  2exp-iq(xi-xi ) ] 

[1  ♦  Z  k,  i  I  b(k,l)  |  2expiq(xa-xi )  ]  +  [1*--*  2) 

-  [1  +  Z  k.i  |  d(k,l)  |  *] [1  +  Z  k.i  |  b(k,l)  |  2]  +  cos(K'-K")(xi-x*) 

+  Zk.il  d(k,l)|  *cos[(K'  -K")(x*-xi)  +  q(xi-xa )  J 

♦  Z  k, i |  b(k,l)  |  2cos[  (K’  -K" ) (xa-xi )  -q(xi-xa)] 

+  Z  kk>  1 1  •  |  d(k,l)  |  >|  b(kM’)|  2cos(K’  -Ku+q,-q)  (xi-xi) .  (7.3) 

Taking  into  account  that  the  momentum  changes  caused  by  infraquanta  are  very 
small  (q<<Ki-K»),  we  now  smooth  this  rapidly  oscillating  correlation  by 
averaging  over  short  distances  of  the  order  l/(Ki-Ki)  and  get  zero,  except 
for  Ki  =  Kj,  which  still  allows  each  of  the  incoming  directions  and  outgoing 
directions  to  be  different 

A(Ki,Ka,K’,K";xi,xa)  =  <S|  0|  S>  ?  [1+Z  k.i  |  /9(k,l)  |  2][1+Z  k.i  |  b(k,l)  |  2] 

+  (1  +  Z  k .  i (  |  /9(k,  1) )  |  2  +  |  b(k,l)  |  2]cos(q(xj-xa  Sjdr  ,«» 


+  £  kh'.ia  ’  |  l(k,l)  I  *|b(kM')|  2cos(q-q’ )  (xt-xi)Sy  ,  i*. 


(7.4) 


This  result  indicates  that  only  particles  of  the  same  energy  have  quantum 
1/f  noise  cross-correlations,  and  groups  of  different  energy  yield  indepen¬ 
dent  quantum  1/f  noise,  as  was  anticipated  empirically  by  Kleinpenning 1 0 . 
Our  result  in  Eq.  (7.4)  differs  somewhat  from  a  heuristic  generalization  of 
the  basic  quantum  1/f  noise  formula  used  in  a  previous  paper*. 

For  fermions  we  start  from  Eq.  (3.7)  and  obtain 


A(Ki,K,,K’,K";xi,x2)  =  ( 1/2) [  |  x(xi)|  *|  tf(x,)|  *  +  (1--2)] 
-(l/4)[**(xi)tf>*(xi);E(xi)<£{xi)  +  <1«--2)1 

3  [1  +Z  k,  i  |  >9(k,l)  |  *]  [1  +  £  k.  i  |  b(k,l)  |  *]  -  ( l/2)cos[  (K'-K* )  (xi-xj) ] 


-(1/2)2  k,.[  I  /9(k,l) |  *  +  |  b(k,l)  |  *)cos[(K'-K")(X2-Xi)  +q(x»-x2)] 
-(1/2)2  kkMfl  /8<k,l)|  *|b(kM')|  *cos(K*  -K"+q'-q)  (xi-x* ) . 


Applying  again  the  smoothing  operation  we  obtain  now 


A(Ki ,K2 ,K’  ,K" ;xi ,x2 )  =  [1+2  k, i  |  jfl(k,l) |  1 ] [1+2  k .  i  |  b(k,l) |  *] 
-<1/2){1  +  2  k .  l  [  |  /9(k,l) )  |  »  +  |  b(k,l)  |  »]cos[q(xi-X2)J>dr.r* 

-(1/2)2  kk* ,  1 1  *  I  /9(k,l)  |  *  |  b(k’,l’)  |  *cos(q-q’ )  (xi-x*)^* 


(7.5) 


(7.6) 


For  elastic  scattering  we  have  K‘  =  Kj  and  K"  and  Kg,  provided  the  scattering 
center  is  fixed  or  very  massive  (e.g.  a  crystal).  We  include  the  infrared 
radiative  corrections  this  time. 

Using  Eq.  (4.16),  the  angular  integrations  can  be  performed  in  Eq.  (7.4). 
We  also  apply  the  relations  dk/k  =  df/f  and  k/k0  *  f/f0  and  obtain 


n 

if 

m. 


!»i:V 


m 

te 

m 


I 

m 

mi 


I 


m 


AbtKpKg.K'.  ^’.Xi.xg)  ■  Cl  +  aAj  j  (f/f0)aAldf/f][l  +  aA2  J  ( f /f 0 ) ^df / f ] 

'o  fo 

+  11  +  aJF[Al(f/f0)aAl  +  A2(f/fo)aA23cosq(xi-x2)(df/f)  (7.7) 

*o 

♦  «2Al«2  JF(W0)°Al(df/f)  |F(f7fo)°A2(df'/f')cos(q-q')(xr,,2>t«K.>K.. 

with  Aj  *  2fl^(K'-Ki)^/3irm^c^  and  A2  ■  2ff2(K"-K2)2^3ffm^c^  for  bosons.  The  lower 
frequency  limit  f0  is  the  resolution  limit  given  by  the  reciprocal  duration  of 
the  noise  measurement.  The  upper  limit  F  is  very  high,  close  to  the  energy  of. 
the  scattered  particles,  divided  by  Planck's  constant.  From  Eq.  (7.6)  we  obtain 

Af(K1,K2,K',  ^",X1,X2)  -  Cl  +  oAx  JF(f/f0)°Aldf/f][l  +  aA2  JF(f/f0)aA2df/f] 

fo  f0 

-(1/2)(1  +  ajF[A1(f/f0)aAl  +  A2(f/f0)aA2]C0Sq(xrx2)(df/f)}6K,>r  (7.8) 

f0 


-(a2/2)AxA2  JF(f/f0)aAl(df/f)  iF(f,/f0)aA2(df,/f,)cos(q-q*)(x1-x2))«Kl^M, 
f0  *0 

for  fermions.  The  Infrared  radiative  correction  factors  (f/f0)aA  were 

calculated  elsewhere® and  can  be  set  »  1  for  all  practical  cases.  Neglecting  the 

"noise  of  noise"  term  proportional  to  a2,  and  using  the  Wiener-Khintchlne 

theorem,  we  obtain  the  corresponding  cross-spectral  densities  of  the  scattering 

rate  fluctuations  aw+  ♦, >  leading  from  Kj  to  K* ,  and  aw*  ♦  >  from  K2  to  K" 

Kj,K  *^2 ,  k 

SM(Kl.K2XX')b  -  (a/2f)[A1(f/f0)aAl+A2(f/f0)aA2]<w^^|><w^^((>5K,tKl,  (7.9) 

for  bosons,  the. factor  in  front  being  a/%<f  in.  general  (N»2  here),  and 


S„(Hi.lCj,K,.K-)f  -  (^f)CA1{f/f0)oAltA2(f/f0)a*2]<>«.  .,><«.  (7.10) 

n2,n  n 

for  fermions,  the  factor  in  front  being  a/(N^>-i)f  now  (see  Sec.  VIII). 


VIII.  PAIR-CORRELATION  FUNCTION  FOR  N  OUTGOING  PARTICLES 

So  far  we  have  considered  only  the  simplest  case  of  two 
particles  present  in  the  outgoing  state,  which  is  the  minimal 
number  that  allows  for  correlations  between  different  particles 
to  be  defined.  Quantum  1/f  noise  is  the  manifestation  of  long- 
range  correlations  between  particles,  causing  them  to  be  bunched 
and  to  be  subject  to  superpoissonian  statistics.  Therefore,  we 
are  interested  in  these  correlations  which  reveal  themselves  in 
the  pair-correlation  function,  also  known  as  the  two-particle 
distribution  function. 

In  general,  there  will  be  more  than  two  particles  in  the 
outgoing  state,  which  are  observed  in  order  to  determine  the 
scattered  current.  In  calculating  the  pair-correlation  function, 
we  should  therefore  consider  the  general  case  in  which  N  scat¬ 
tered  particles  are  simultaneously  under  observation  in  the  out¬ 
going  state  during  the  measurement  of  the  cross  section.  This 
is  what  we  will  try  to  do  in  tis  section. 

Similar  to  the  case  of  two  particles  in  Eq.  (2.1).  the 
state  of  N  scattered  particles  can  be  written  in  the  form 


Here  we  have  denoted  by  ?i(£i)  the  exact  single-particle  wave  function  of 
the  particle  i,  where  both  the  interaction  with  the  scattering  center  and 

mtiih  j 

the  interaction 'the  electromagnetic  field  modes  are  taken  into  account.  For 
bosons  the  operator  of  the  pair-correlation  is  given  by  Eq.  (3.3),  and  by 
using  the  permutation  relations  Eqs.  (3.2)-(3.2b),  we  obtain 


N!  <S°  |  0|  S»>  =  <0J  ^(^N)...^(>?i)0+(Xi)^*(xj)|*{xa)^(xi)^(5i)...^(fii)  |  0> 
=  <0|in  iHVi)d(VY  -Xi)S{t?a-X2)Z  5  \  o> 


=  2  $(^-Xi)d<77/i-xa)<S(£B-Xi)$(f.-x*)2  .  Il'  Ij'  S(ViSi) 

=  22  $(9*-xi)$(»A-x*)£(£b-xi)£(£.-x*)2  II  diVi-Si), 


(8.2) 


where  the  sun  over  permutations  runs  over  all  permutations  of  the  N-2  in¬ 
dices  j.  This  allows  us  to  calculate  the  complete  matrix  element 

f/  ff 

<S  |  O  |  S>  =  [1/N(N-1)]2  2  Jd’nA  Jd*» Td»£«Jd3£B 

fiM  «***»•!  T  '  (y 

^*(^)^  Mjl^)^(f*)<M £■)<?( *7^ -XiW^-XjJS^i-XjJd^.-Xj) 

=  [1/N(N-1)]2  2  th.  *(xa)^  «(xi)tf.(xi)tfB(x*) 

\  /  \  N 

=  [1  ♦  2  k.i  |  b(k,  1)  |  *]-MN(N-l)  +  [ (N-l)/N) [22  -j2  k.i  I  b.(k,l)  |  * 
aJ 

+  22  «<i  2  k.i  |  b«(k,l)  |  *cosq(xi-xj)]J 


=  tl  +  2  k.i  1  b(k,l)  J  *]-MN(N-l)  +  2(N-1) [2  k.i  I  b(k,l)  |  * 
+  2  k, i  |  b(k,l)  I  *cosq(xt-Xj)3) . 


(8.3) 


Dividing  the  term  dependent  on  xi-xj  by  the  constant  term,  and  neglecting 
terms  of  higher  order  than  |  b(k,l)  |  * ,  we  obtain  the  fractional  spectrum 
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density  of  the  quantum  1/f  fluctuations  in  the  concentration  e  of  particles, 
in  the  current  density  j,  and  in  the  physical  cross  section  <3 


Sj(f)/<j>*  =  2oA/Nf. 


(8.4) 


In  Eq.  (8.3)  we  have  used  the  form  of  the  single-particle  wave  functions 
given  by  Eq.  (5.1)  or  (5.2),  and  in  Eq.  (8.4)  we  have  used  the  expression  of 
the  bremsstrahlung  amplitude  given  by  Eq.  (4.16).  Our  calculation  in  Eq. 
(8.3)  and  the  result  in  Eq.  (8.4)  show  how  the  independent  sets  of  phases  in 
the  bremsstrahlung  energy  loss  part  of  the  single-particle  wave  functions 
lead  to  independent  additive  quantum  1/f  spectral  contributions  from  each  of 
the  scattered  particles,  and  to  the  1/N  factor  in  the  fractional  spectral 
density  shown  in  Eq.  (8.4).  This  result  and  the  foregoing  calculation  are 
in  agreement  with  our  previous  results1**  and  calculations,  and  perfectly 
justify  our  previous  use  of  single-particle  wave  functions  with  subsequent 
introduction  of  the  1/N  factor  on  the  basis  of  the  independence  of  1/f  noise 
contributions  from  different  carriers. 

In  the  case  of  fermions  the  calculation  is  similar,  except  for  the 
use  of  anticcmnutators  for  the  fermion  field  operators.  In  order  to  em¬ 
phasize  the  independence  of  our  results  on  the  representation  used,  and  to 
show  directly  how  the  calculations  presented  in  this  paper  can  be  performed 
without  second  quantization,  we  give  here  the  direct  calculation  in  terms  of 
a  Slater  determinant  for  the  state  of  N  scattered  fermions 


tfi i , . . . m(ri. . .r»)  =  (N!)*1'1  I. 
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Here  r  combines  the  position  vector  x  and  the  spin  variable  s  for  each  of 
the  particles.  The  pair-correlation  function  is  obtained  by  integrating 
with  respect  to  the  coordinates  of  all  but  two  of  the  fermions 

A(ri,r2)  -  Jd*r3« .  .d3r»<  |  1 . . .  t  «(ri. .  .m)  |  J>«  (8.6) 

Here  the  integrals  also  include  summations  over  the  spins,  and  the  expecta¬ 
tion  value  is  with  respect  to  the  phases  present  in  the  bremsstrahlung  parts 
of  the  wave  functions  and,  if  the  emitted  photons  are  included  in  the  final 
state,  the  expectation  value  is  also  done  over  the  vacuum  of  the  photons. 
Assuming  orthonormality  of  the  functions  4> n...<£in,  we  obtain 

H 

A(r»,r»)  =  [1/N(N-1)]2  a,..i<  |  -  *.(r»)*,(ri)  |  *>.  (8.7) 

To  display  the  spin  variables  explicitly,  we  write  ^a(r)  =  Xm(x)  |  s>  and  get 

H*. 

A«t'(xi,xa)  =  [1/N(N-1)]Z  ■,b«i<[x«*(xiKS|  Z»*(x*)<s’ |  (8.8) 

-  Xmt{Xi)<B>  I  Xb*(Xi)<S  I  ]  [*.(Xi)|s>  X,(xa)|s’>  -  X«(X2)|s'>  Z.(xi)fs>  ]> 

HA 

=  [1/N(N-1)]Z  |  Xa(x»)  |  1 1  ZB(X2)  |  *  +  |  x.(x2)  I  *|  *.(x0  |  * 

-  *.*<Xi)*.(Xj)  I  <s|  s’>|  3Xb*(x* )Xb(Xi)  -  (Xi«--*X2)]>. 

Here  the  simbol  (xi*-  -*•  x8)  designates  the  immediately  proceeding  term.  Con¬ 
sidering  all  spin  orientations,  we  obtain 


A(xi ,xi)  - 


H/2 

A  ♦  A  =  (1/  N(N-1 ) ] 2  ■,b.i<(4  |  Xm\  3  I  |  3 

n  r* 


& 


-  Xra*(xi)x*n(x2)xm(x2)xn(xl)  -  (xx  ♦  ♦  x2)]>  (form  1) 

»  C|C|  4/x2N(N-l)3lN2[l+EjT i|b(k,l)|232 

-2  E  exp[1Kffl(xi-X2)3Cl  +  2jc,l  lbm(k»1 )  |z  expClqm(xi-X2)3 
m*l 
N/2 

•E  exp[-1Kn(x1-x2)3Cl  +  •  i  •  |bn(k' ,1  *)  I2  exp[1qn(x1-x2)3  (form  2) 

n*l 

2  [|C|  4/x2(N-l)3iNCl  +  E|c,iib(k,l)|  232 

N/2  -  _ 

-  (2/N)  E  [1  +  2£ff  t|bn(k,l)|  2cosqn(x1-x2) 

n*l 

+*  |bn(lc,l )  |2  |bn(iT'  ,1 ')  |  2cos(qn-q,n)(x1-x2)3} 

k.l;k\r 

*  C jc I  4/x2(N-1)3(N[1  +  Ej?j|b(M)|  z32  -1  -2EjT>1|b(iT,l)|  ^osq^-Xg)}.  (8.9) 

Hiis  form  of  the  pair-correlation  function  includes  the  1/N  factor  which 
multiplies  the  variable  (noise)  part.  The  crucial  point  in  the  derivation 
of  the  1/N  factor  was  a  elimination  of  the  rapidly  oscillating  terms  exp(KB- 
K*)(xi~xt)  with  Ka  5*  Ka  present  in  the  second  form  of  Gq.  (8.9)  above,  an 
elimination  indicated  through  the  approximation  sign  connecting  the  second 
form  to  the  third  form  above.  Indeed,  since  Ka  differs  from  K*  by  much  more 
then  the  momentum  change  corresponding  to  the  emission  of  a  infraquantun, 
these  terms  will  have  a  very  fast  oscillation,  and  will  not  yield  any  low 
frequency  noise.  Since  they  are  also  small  in  magnitude,  they  are  negli¬ 
gible.  This  provides  the  important  reduction  of  the  noise  term  by  a  factor 
N.  In  the  last  form  of  Eq.  (8.8)  the  first  two  terms  are  constant  and 
large,  and  do  not  yield  any  rapid  oscillations  which  would  justify  elimina¬ 
tion  of  any  cross  terms. 

In  Eq.  (8.9)  we  have  used  again  the  form  of  the  single-particle  wave 
functions  given  in  Eq.  (5.1)  or  (5.2),  with  independent  sets  of  phases 
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present  in  the  bremsstrahlung  energy  loss  parts  of  each  particle.  Had  we 
used  a  random  time-shift  in  each  of  the  single-particle  wave  functions,  i.e. 
just  a  random  initial  time  constant,  or,  equivalently,  a  random  space-shift- 
constant,  all  results  would  have  been  exactly  the  same.  In  fact,  the  random 
shift  should  better  describe  the  initial  Poisson  distribution  of  the  incom¬ 
ing  particles  which  are  scattered.  As  we  have  seen  in  Sec.  IV,  the  random 
phase  set  in  the  bremsstrahlung  energy  loss  parts  is  the  set  of  random  ini¬ 
tial  phases  of  the  electromagnetic  field  oscillators.  Therefore  it  should 
come  in  the  same  way  for  all  particles.  However,  the  random  shift  will 
eliminate  expectations  of  the  cross  products  b»*(k,l)ba(k,l)  with  m  *  n 
just  as  the  sets  of  random  phases  used  by  us  did;  if  p  is  for  instance  a 
random  space  shift,  these  cross  terms  yield  contributions  to  the  pair- 
correlation  function  of  the  form  <|  b(k,l)  |  *  cosq(xi-Xj  +  p)>  =  0,  where 
the  average  is  with  respect  to  p.  We  present  this  observation  here  as  an  af¬ 
terthought,  because  the  sets  of  random  phases  generated  by  a  shift  in  time 
or  space  will  appear  to  be  random,  but  will  still  contain  some  correlations. 
The  point  we  make  here  is  that  these  correlations  have  no  effect  on  our  cal¬ 
culations,  so  we  can  continue  to  use  random  phases  for  our  purpose.  There 
may  be,  of  course,  some  differences  in  higher-order  correlations  which  we  do 
not  consider  here. 

In  the  last  form  of  Eq.  (8.9)  we  have  neglected  the  higher-order 
term.  Using  Eq.  (4.16),  we  can  write  the  pair-correlation  function  for  fer¬ 
mions  in  the  form 


A(x!,x2)  -  (|C|  4/x2)  11  +  [2/(N-l)]ieil|b(ic,l)|  2[Nct,  cosq(xrx2)]}.  (8.10) 
uividing  again  the  variable  part  by  the  constant  terra,  and  neglecting  small 
constant  terms,  we  obtain  for  the  fractional  spectral  density  of  the  fermion 


current  and  cross  section  fluctuations 


Sj(k)dk/j2  -  2oAdk/k(N-l)  -  Sj(f)df/j2  -  2oAdf/f (N-l).  (8.11) 


IX.  DISCUSSION 


As  we  have  seen  in  Eqs.  (2.5a)  and  (2.5b)  the  second-quantized  cal¬ 
culation  of  the  pair-correlation  function  for  bosons  has  yielded  just  the 
autocorrelation  of  the  single- particle  probability  density  which  we  had  used 
in  the  initial  publications 1 • * .  For  fermions  the  same  result  was  obtained 
with  a  minus  sign  which  corresponds  to  a  180*  phase  shift.  The  latter  was 
interpreted  as  the  result  of  the  Pauli  exclusion  principle  which  prohibits 
the  very  close  location  of  two  identical  fermions  of  the  same  spin. 

However,  this  affects  the  pair-correlation  function  only  at  very  small  rela¬ 
tive  distances  of  the  particles.  For  the  larger  distances  and  lower 
frequencies  or  wave  nunbers  important  for  1/f  noise,  this  difference,  which 
appears  in  the  case  of  fermions,  is  negligible.  We  conclude  that  the  simple 
autocorrelation  function  calculated  for  the  single-particle  probability  den¬ 
sity  in  earlier  papers  on  quantum  1/f  noise  is  required  by  any  second- 
quantized  calculation  which  takes  into  account  the  quantun  exchange  between 
identical  particles  automatically.  However,  there  is  nothing  special  about 
the  second  quantization  beyond  this  automatism  and  its  neatness;  the  same 
results  are  obtained  also  without  it,  by  properly  symnetrizing  the  wave 

functions,  as  we  have  seen  directly  in  Sec.  VTII..  Furthrrmore,  quintum  i/f 
noise  m  also  derived  without  Che  APSPSWF. 
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At  first  sight  the  autocorrelation  of  the  single-particle  probability 


density  function  can  not  be  considered  as  the  expectation  value  of  an  ob¬ 


servable  quantity  in  agreement  with  the  principles  of  quantum  mechanics,  be¬ 


cause  it  is  of  fourth  order  with  respect  to  the  wave  function.  However,  as 


we  have  seen,  this  is  required  by  a  correct  calculation  in  second  quatiza- 


tion.  To  explain  the  fourth -order  dependence  of  the  pair-correlation  func¬ 


tion  on  the  wave  function  without  using  the  method  of  second  quantization, 


we  mention  that  the  nonuniformity  of  the  distribution  of  the  outgoing  par¬ 


ticles  is  described  by  the  two-particle  wave  function,  which  in  turn  is  a 


properly  symmetrized  or  antisymmetrized  product  of  two  single-particle  wave 


functions.  Taking  the  module  squared  of  this  two- particle  wave  function^ we 
obtain  the  probability  density  of  the  relative  positions  of  the  par¬ 


ticles,  or  the  pair-correlation  function,  which  therefore  will  be  of  fourth- 


order  in  terms  of  the  single-particle  wave  functions.  This  fourth-order  de¬ 


pendence  which  was  present  in  all  publications  on  quantum  1/f  noise1”*,  is 


therefore  in  agreement  with  the  principles  of  quantum  mechanics,  and  is  jus¬ 


tified  by  the  detailed  many-particle  calculation,  in  spite  of  the 


acrimonious  criticism  it  has  recieved.  We  have  also  provided  for  the  first 


time  in  this  paper  a  formal  derivation  both  of  the  quantum  1/f  cross  cor¬ 


relations  and  of  the  1/N  factor  of  the  empirical  Hooge  relations,  in  perfect 


agreement  with  the  previously  introduced  1/N  factor1-*. 


In  conclusion,  we  have  shown  how  the  method  of  second  quantization 


can  be  applied  to  derive  the  pair-correlation  function  of  the  conventional 


quantum  1/f  effect  in  space  and  in  time.  Both  the  pair-correlation  function 


and  the  power  spectral  density  are  similar  to  the  previously  published  ex¬ 


pressions  based  on  more  elementary  derivations.  Rather  than  contradicting 


the  existence  of  quantum  1/f  noise,  the  second  quantization  method  confirms 


the  quantum  1/f  effect  in  a  brilliant  way,  and  convincingly  justifies  my 
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previous  use  of  four  single-particle  wave  functions  in  order  to  describe  the 
self-interference  of  the  particle  fields,  on  the  basis  of  quantum  exchange 
effects  between  identical  particles. 

Finally,  we  mention  that  the  inclusion  of  the  soft  photons  into  the 
final  state  is  not  required,  if  the  electromagnetic  field  inodes  are  not  part 
of  the  system  studied,  but  are  considered  external,  being  included  as  a 
time-dependent  external  force  field  in  the  Hamiltonian.  This  external  field 
description7  •  •  involves  a  time-dependent  Hamiltonian  and  a  formal  lack  of 
energy  conservation,  but  is  well  suited  for  the  derivation  of  quantum  1/f 
noise  and  of  most  other  electrophysical  and  electronics  problems,  as  I  have 
pointed  out  earlier.  The  present  paper  shows  that  the  quantum  1/f  effect 
can  also  be  derived  in  the  language  of  second  quantization,  with  the  ex¬ 
plicit  inclusion  of  the  emitted  bremsstrahlung  photons  into  the  final  state. 

The  equivalence  of  the  wave  number  spectrum  with  the  frequency 
spectrum  was  proven  here  by  neglecting  coherent  state  quantum  1/f  noise, 
which  was  discussed  elsewhere  * *• *  * .  Indeed,  if  the  coherent  state  quantum 
1/f  noise  is  included,  fluctuations  caused  by  virtual  photons  in  the  final 
state  must  be  included.  We  recall  that  coherent  quantum  1/f  noise  is  caused 
by  the  energy  uncertainty  introduced  by  the  coherent  state  of  the 
electromagnetic  field  of  a  physical  charged  particle.  Both  coherent  and  in¬ 
coherent  quantum  1/f  noise  can  be  derived  in  the  interaction  picture.  In 
the  interaction  picture  coherent  quantum  1/f  noise  is  obtained  from  diagrams 
which  do  not  involve  scattering  of  the  particles,  while  the  conventional  (or 
incoherent)  quantum  1/f  noise  considered  in  this  paper  arises  from  diagrams 
of  scattering  with  the  simultaneous  emission  of  photons.  We  have  not  in¬ 
cluded  coherent  quantum  1/f  noise  in  our  calculation  here,  because  we  have 
considered  only  the  noise  in  the  final  state  long  after  the  scattering 
process,  with  the  coupling  to  the  electromagnetic  field  disconnected,  and 
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with  the  Heisenberg  field  operators  <A( r,t)  in  free  motion.  Without  this 
restriction,  the  quantum  1/f  noise  spectrum  in  frequencies  would  have  been 
larger  than  the  wave  number  spectrum.  However,  in  very  small  electronic 
devices  only  conventional  quantum  1/f  noise  seems  to  be  present M. 
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APPENDIX  ft 


THE  SINGLE-PARTICLE  WAVE  FUNCTIONS 

1.  Introduction 

The  single-particle  wavefunctions  used  in  this  paper  make  sense  physically 
and  were  first  Introduced  in  1975  on  physical  grounds*.  The  objective  of  this 
appendix  is  to  show  how  they  can  be  derived  formally.  The  single-particle  wave 
functions  were  derived  in  Sec.  IV  by  treating  the  electromagnetic  field 
classically,  and  by  writing  the  final  result  In  terms  of  photon  creation  and 
annihilation  operators.  This  last  step  implies  the  correspondence  principle. 

A  more  rigorous  derivation  should  follow  perturbation  theory,  e.g.  the 
calculation  of  Yennie,  Frautschi  and  Suura*®  reformulated  by  Chung*®  In  order  to 
provide  a  scattering  matrix  element  free  of  infrared  divergences.  The 
reformulation  by  Chung  shows  that  the  price  we  have  to  pay  for  a  finite  matrix 
element  Is  a  more  complicated  set  of  Incoming  and  outgoing  states  describing  the 
scattered  particles,  with  the  inclusion  of  a  coherent  state  of  the 
electromagnetic  field  Into  the  notion  of  charged  particle.  This  new  picture  was 
first  suggested  by  Dollard*7  In  the  language  of  nonrelatl vlstlc  quantum 
mechanics,  and  was  successfully  generalized  to  the  relativistic  case  by  Kulish 
and  Faddeev*^,  and  by  Zwanzlger*^.  This  new  picture  is  intermediary,  between 
the  usual  Heisenberg  picture  and  the  usual  interaction  picture,  because  the 
offending  (long-range)  part  of  the  electromagnetic  interaction  has  now  been 
included  into  the  unperturbed  hamiltonian. 

Following  Chung*®,  we  obtain  a  finite  matrix  element  by  writing  the  initial 
state  of  a  charged  particle  involved  in  a  scattering  process,  with  initial 
4-momentum  p^  and  final  4-momentum  pf,  in  the  form 

|  >1  -  expt-(l/2)rktJl|S-^(k)  |2} 

exPlSk.i  s{*\k)e(*)(k)all(A)+)  |^(Pi)> 


•  U  -  (l/2)£lt,t|si‘)(k)|2nl  +  Ik,t  S1U)(k)e(*)(k)«|tU)+1  |V(P1)>,  (Al.l) 
where  the  summations  over  k  can  be  transformed  In  Intergrals  by  multiplication 
with  (L/2w)3  as  In  Eq.  (5.4),  and  where 

SiU)(k)  -  (e/2(2it)3k0)p1e(^)/kpi  (A1.2) 

Is  a  function  which  depends  on  the  momentum  pj  of  the  Incoming  particle  of 
charge  e.  The  state  | ( p-| ) >  includes  the  bare  particle  and  the  vacuum  on  which 

(l)+  fa\ 

the  creation  operates  a^  act.  Here  e^*'(k)  Is  the  polarization  4-vector  of 
the  mode  k,i.  The  physical  charged  particle  defined  by  the  first  form  of  Eq. 
(Al.l)  contains  the  bare  particle  and  Its  field  which  Is  in  a  coherent  state, 
with  an  Indefinite  number  of  photons  In  each  mode  and  therefore  with  energy  and 
momentum  which  are  not  sharp,  l.e.,  not  well  defined.  The  last  form  of  (Eq. 
Al.l)  is  only  an  approximation  obtained  by  expanding  the  exponentials.  The  *zrm 
of  5j(*)(k)  given  by  Eq.  (A1.2)  is  correct  only  in  the  limit  k+0  of  Interest  for 
us. 

The  final  state  is  obtained  by  replacing  p^  with  pf  In  Eqs.  (Al.l)  and 
(A1.2).  The  photon  Is  assumed  to  have  a  finite  mass  \  which  Is  allowed  to 
approach  zero  at  the  end  of  the  calculation. 

It  Is  easy  to  verify  that  if  the  initial  and  final  particle  states  are 
chosen  as  shown  In  Eqs.  (Al.l)  and  (A1.2),  all  infrared  divergences  cancel 
already  at  the  level  of  the  matrix  element  Mf ^ ,  and  not  only  later  at  the  level 
of  the  process  rate  or  cross  section  as  would  be  the  case  had  we  started  with 
the  usual  notion  of  particle,  without  the  long-range  part  of  the  field  included. 
It  is  easier  to  do  this  proof  in  two  steps:  first  to  second  order  only,  in 
order  to  get  acquainted  with  the  formalism;  then  up  to  any  order.  Our  treatment 
will  be  simplified  and  schematized  as  much  as  possible,  in  order  to  focus  on  the 
essential . 
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2.  Cancellation  of  Infrared  divergences 


Infrared  divergences  are  of  two  kinds,  arising  from  virtual  photons  and  real 
photons.  The  former  are  photons  which  boomerang  back  to  one  of  the  charged 
particles  Involved  in  the  contemplated  process.  Real  photons  are  emitted  (or 
absorbed)  for  good.  Let  Mq  be  the  basic  matrix  element  calculated  for  an 
arbitrary  process  (e.g.,  scattering  of  charged  particles  on  a  fixed  potential) 
without  Including  the  coupling  to  photons.  In  fact,  however,  the  process  may 
also  happen  with  any  number  of  virtual  photons.  Adding  up  the  Infinite  series 
of  processes  (or  Feynman  diagrams)  corresponding  to  virtual  processes  of  any 
order,  one  obtains^**®  the  matrix  element 


M'o  »  V018  . 

where  we  have  neglected  a  small  nondlvergent  term,  and 


(A1.3) 


Re(cxB)  =*  (e2/4V)z(k2+X2)1/2  - - 

k  U2Pf 


(2(pf-k)  (2(p,-k)-](2) 

k-x?)  -  (2p,.k-x)J  <A1-4> 


-  -a/2)ik.tlSf‘,(k)-SlU,(MI2  . 

Here  the  rectangular  bracket  Is  a  4-vector  and  the  exponent  In  brackets 
designates  the  scalar  product  of  the  four-vector  with  Itself.  For  x  »  0  the 
Integral  derived  from  the  sum  in  Eq.  (A1.4)  shows  a  logarithmic  divergence  which 
has  been  exponentiated  In  Eq.  (A1.3)  by  summing  the  the  virtual  photon  processes 
(diagrams)  of  all  orders.  The  value  -«  of  Eq.  (A1.4)  tells  us  that  M'  *  0  in 
the  limit  X  *  0,  and  therefore  there  are  no  processes  without  the  simultaneous 
emission  or  absorption  of  any  (observable  or  unobservable)  real  photons. 

To  get  a  finite  result,  we  multiply  the  matrix  element  M0  which  includes  all 
virtual  photon  contributions  with  the  sum  of  all  real  photon  contributions.  To 
second  order,  this  yields 
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M*  -  Hoe08  1 1— ( 1/2)  £  |s{l)(lc)l2Hl-(l/2)  £  IS^’tk)!2] 


11  ♦  £  s}*)tk)s}*)(k)nl  +  £  sjl)(k)su)(k)j^l  -  £  s{t)(k)s(iJk) j  ♦  nMt 

k  f  t»  k  ft 


■  M0[l  +  n  ♦  aB  -  (1/2)  E(|Si  |2  +  |Sf  |2  -2S^Sf)  +  £  |Sk  |2| 
■MoLl+n+aB+  oBj. 


(A1.5) 


Here  In  the  first  form  on  the  r.h.s.  we  have  first  Included  the  normalization 
factors  present  In  the  In  and  out  coherent  states  (Al.l).  The  following  factor 
(third  rectangular  bracket)  Is  the  scalar  product  of  the  In  and  out  coherent 
fields,  describing  a  transition  without  Interaction.  The  next  factor  describes 

U) 

the  emission  of  the  final  coherent  state  field  Sf  (k)  In  the  scattering 
process,  and  the  last  factor  corresponds  to  the  absorption  of  the  Incoming  state 
coherent  field.  Note  that  the  last  three  rectangular  brackets  have  just  been 
taken  over  and  copied  from  the  usual  treatment  of  infrared  divergences,  just  as 
we  took  over  the  virtual  photon  contribution.  In  the  second  form  on  the  r.h.s. 
we  consequently  limited  ourselves  to  the  second  order  In  e/(hc)^2  and  carried 
out  the  multiplications.  We  omitted  the  upper  polarlzatln  label  and  wrote  the 
argument  k  as  a  superscript,  also  omitting  the  Indication  of  the  summation 
variables  k  and  i  under  the  sign  of  the  sum.  Finally  we  Introduced  the 
notations 

Sf4)(k)  -  SiU)(k)  =  sU)(k)  =  S*  ;  aS  -  (l/2)z|S^|2  ,  (A1.6) 

In  order  to  save  space.  The  term  M0n  stands  for  a  non-inf rared-di vergent 
contribution.  We  notice  from  (A1.4)-(A1.6)  that  for  low  k  values  the  summands 
(or  Integrands)  In  aB  and  aB  are  Identical,  but  with  opposite  signs  in  the  limit 
A-*0.  Therefore,  the  infrared  divergences  cancel  already  at  the  level  of  the 
matrix  element. 
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The  coherent  states  of  the  field  of  a  physical  charged  particle  defined  by 
Eq.  (Al.l)  are  vectors  In  a  separable  Hilbert  space.  The  result  In  Eq.  (A1.5) 


Is  not  affected  by  a  translation  In  this  space,  which  amounts  to  adding  or 

subtracting  the  same  vector  yk  from  all  states.  Indeed,  subtracting  y*f  the 
Ic  Ic 

amplitude  Sk  *  Sf  -  Sj  remains  the  same,  and  we  get 

M“fi  a  MoeaBl1-(1/2)s|Si-yk|2jLi-(i/2)E|Sf-yk|2j 

k  k  k  k 

Ll-ReE(Si-yk)*(Sryk)j[l+Rez(Sf-yk)SkJtl-Rez(S1-yk)Skj  +  nM0  . 


a  M0(l+n+aB+aB) 


(AI.7) 


as  before. 


To  demonstrate  the  cancellation  of  divergences  to  all  orders,  we  regroup  the 
perturbation  theory  series  in  the  form 


k  +  + 

Mf1  *  Mo<0|expLSfakjexpi>zSk(ak-ak)Jexp[-£SialcJ  |0> 

expt-(l/2)r(|S?|2  +  |Sfl2)j  ♦  Hon' 
where  Mnn'  is  a  nondlvergent  part. 


(A1.8) 


Applying  the  relation  e^+^  *  ePeQe-(l/2)[P,Q]  with  P  a  zSkak  and  Q  *  P+,  we 
obtain 

k  k  kb 

Mf i  »  MoeaB<0|expcSfakexpESiak|0><0|expESfakexpESkak|0> 

<0  lexp-CS^explS^IO  *  K„n'  ■  ♦  „•).  (A1.9) 

The  virtual  phonon  factor  e°®  arises  from  e"^2^P»^]  automatically,  and  has  an 
upper  Integration  limit  roughly  equal  to  the  energy  E  of  the  particle  considered 
and  involved  in  the  scattering  process.  The  integrals  in  eoB,  however,  should 
be  divided  in  a  part  below  the  photon  observation  threshold  e0  *  hf0  and  a  part 
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above,  going  also  up  to  E,  which  contains  the  part  of  the  matrix  element 
occuring  with  bremsstrahlung.  The  parts  below  c0  cancel  each  other  In  aB  and  In 
aB,  demonstrating  the  disappearance  of  all  divergences.  Bringing  the  rest  of 
e°®  to  the  denominator  and  neglecting  the  difference  between  B  and  §,  we  can 


write 


1  -  [expi'IS*!2  -  lj 
^  expz*  |Sk |2 


(A1.10) 


The  sum  with  a  prime  is  over  |k |  >  e0  only,  i.e.,  is  restricted  to  observable 
photons.  Here  we  have  separated  in  the  numerator  the  unity  which  represents  the 
non-bremsstrahlung  part  of  the  rate,  as  we  did  earlier*  to  second  order. 


3.  Single-particle  wave  functions 

The  discussion  and  derivation  outlined  above  suggests  Eq.  (Al.l)  as  the 
correct  choice  for  the  incoming  and  outgoing  physical  charged  particle  states. 

For  the  out  state  the  subscript  i  is  replaced  by  f.  However,  conventional 
quantum  1/f  noise  (effect)  is  defined  as  the  cross  section  (or  process  rate) 
fluctuation  arising  due  to  the  small  bremsstrahlung  effects  associated  with  the 
process  considered.  The  current  fluctuations  caused  by  the  general  uncertainty 
in  the  energy  and  momentum  of  the  physical  charged  particle  states  have  always 
been  called  coherent1***^  state  quantum  1/f  noise  (effect),  and  will  be  studied 
separately  in  another  paper**»*^.  Since  we  are  interested  only  in  the  noise 
contribution  introduced  by  the  collision  at  hand,  we  perform  a  translation  in 
the  space  of  principal  (coherent)  state  vectors,  by  taking  y  *  in  Eq.  (A1.7). 

The  Invariance  of  matrix  elements  with  respect  to  any  translation  was  verified 
above  In  second  order,  and  is  manifest  in  Eqs.  (A1.10)  and  (A1.9)  to  any  order, 

f>  k 

because  S*  is  an  invariant  difference  of  two  vectors.  Subtracting  S-j  from  the  lkl>£® 
amplitudes  of  all  vectors  we  obtain  a  bare  particle  incoming  state  with  no  obsttrvaile 
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photons  (|k|>eQ),  and  an  out  state 


l»f  -  exp{-(l/2)i?>ll|S*(k)|2} 

’  tl-(l/2)^.tIS*(kmitIk-its“)(k)e-,P(I-?,a^)+]|r(Pf»(A1-ii) 


Here  we  have  introduced  also  the  recoil  phase  factors  e"‘'cl(x_r)  which  describe 
the  momentum  loss  q  *  kc/v  of  the  outgoing  paricles  of  momentum  p^  =  mv  when 
a  photon  of  momentum  k  is  emitted.  These  recoil  phase  factors  are  neglected  in 
the  independent  boson  model  which  ignores  the  changes  caused  by  the  emission 
of  individual  photons;  we  have  neglected  these  factors  so  far,  but  are  restoring 
them  here.  The  particle-specific  constant  r  has  to  be  introduced  because  the 
single-particle  states  |>^  are  not  exact  momentum  eigenstates,  and  are  there¬ 
fore  not  exactly  invariant  with  respect  to  translations.  Individual  single 
particle  wave  functions  may  differ  by  such  a  random  translation  in  space  or 
time  regardless  of  whether  they  are  localized  or  completely  unlocalized  in 
space  and  time.  For  momentum  and  energy  eigenfunctions,  such  translations 

introduce  only  a  general  phase  factor  which  drops  out  of  any  physically  relevant 

In  Ej,.  (Al.tt)JTJoitl'  n,r  cirep  cut.j  j  ibCji  r  inf/if  radial  cocr«//«»afc  x  in  ?Ar  n*roneaT,  cgui nleA ft, 

calculationVFor  a  spherical  wave  the  translation  is  replaced  byV^  time-shift, 

or  time  delay;  if  this  requires  more  than  a  general  phase  factor,  a  random ndi a l(cr 

time)shift  must  be  introduced  for  each  scattered  particle  in  the  etpone  nts, 

With  the  notation  S^(k)e^Vqr  =  b(£,t),  Eq.(Al.ll)  becomes  similar  to 

Eq . ( 5 . 2 )  which  was  used  in  the  main  text  of  this  paper.  The  main  difference 

is  the  normalization  factor  in  (Al.ll),  which  was  suppressed  in  Eq.(5.2).  The 

phase  factor  eiqr  (or  elqr  for  spherical  symmetry),  with  constant,  particle- 

specific  r,  present  in  b(l<,t),  is  the  random  phase  factor  introduced  by 


us  in  the  main  text  for  each  particle,  which  led  to  the  1/N  factor  in 


Sec.  VIII.  This  derivation  of  the  single-particle  wave  functions 
treats  the  electromagnetic  field  quantum-mechanical ly .  while  the  deri¬ 
vation  in  Sec.  IV  was  semiclassical . 

In  the  present  paper  we  have  considered  only  one  isolated  scatte- 

L, 

ring  event.  We  also  note  that  the  subtraction  of  Sj  changes  the  noise 
given  by  the  final  state.  The  corresponding  separation  in  conventio¬ 
nal  (Incoherent)  and  coherent  state  quantum  1/f  noise  seems  to  be  con 
venient  both  for  conceptual  or  didactical,  and  for  practical  or  calcu- 
lational  reasons,  connected  with  the  method  of  compounding  the  noise 
contributions  from  successive  scattering  events.  As  was  indicated 
earlier*0,  this  addition  of  contributions  from  many  successive,  closely 
spaced,  scattering  events  gives  a  result  close  to  the  contribution, 
calculated  here,  from  a  single  representative  scattering  event  with  a 
correction  factor  which  is  usually  close  to  unity. 
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PART  II 


APPLICATION  OF  THE  QUANTUM  l/f  THEORY  TO  MIS  INFRARED  DETECTORS 
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I.  INTRODUCTION 


The  utilisation  of  infrared  detectors  at  larger  wavelengths  and  in  the  staring 
mode  has  made  the  limitations  imposed  by  1/f  noise  on  the  performance  of  MIS 
structures  more  stringent  and  conspicuous.  On  the  other  hand*  the  develoDment 
of  the  quantum  1/f  theory  has  provided  us  for  the  first  time  with  a  possibility  to 
predict*  model  and  calculate  from  first  principles  the  various  1/f  noise 
contributions  affecting  the  many  currents  and  processes  which  are  of  importance 
in  the  operation  of  infrared  detectors. 

The  mam  purpose  of  the  present  report  is  to  apply  this  new  knowledge  of  1/f 
noise  to  MIS  structures  working  as  infrared  detectors.  We  start  here  with  a  brief 
general  description  of  quantum  1/f  noise  and  continue  in  Sec.  II  with  the 
inventory  of  various  components  present  in  the  current  Df  MIS  infrared  detectors. 
In  Sec.  Ill  we  analyse  the  quantum  1/f  noise  associated  with  the  currents  and 
processes  discussed  in  Sec.  II.  In  Sec.  IV  we  compare  the  i/f  noise  components  in 
magnitude  and  determine  their  impact  on  the  performance  of  MIS  detectors. 
Finally*  in  Sec.  V  we  discuss  the  resulting  i/f  noise  limitations*  and  point  out 
some  possibilities  of  reducing  the  quantum  1/f  noise.  Those  familiar  both  with 
quantum  1/f  noise  and  with  infrared  detectors  may  now  go  directly  to  Sec.  III. 

Quantum  1/f  noise*-^  is  a  fundamental  fluctuation  of  physical  cross  sections 
and  process  rates*  caused  by  the  infrared  -  divergent  coupling  of  current  carriers 
to  low  frequency  photons  and  other  infraquanta.  The  physical  origin  of  quantum 
1/f  noise  is  easy  to  understand.  Consider  for  example  Coulomb  scattering  of 
electrons  on  a  center  of  force.  The  scattered  electrons  reaching  a  detector  at  a 
given  angle  away  from  the  direction  of  the  incident  beam  are  described  by 
DeBroglie  waves  of  a  frequency  corresponding  to  their  energy.  However,  some  of 
the  elecrons  have  lost  energy  in  the  scattering  process*  due  to  the  emission  of 
Bremsstrahlung.  Therefore*  part  of  the  outgoing  DeBroglie  waves  is  shifted  to 
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slightly  lower  -frequencies  and  inter-feres  with  the  main,  non-Bremsstrahlung, 
component,  yielding  beats.  These  beats  present  in  the  probability  density  along 
the  direction  o-f  the  scattered  beam  will  be  noticed  in  the  detector  as  low 
-frequency  current  fluctuations,  and  will  be  interpreted  as  fundamental  cross 
section  fluctuations.  Although  the  wave  function 'p  of  each  carrier  is  split  into  a 
Bremsstrahlung  part  and  a  non-Eremsstrahlung  part,  no  quantum  i/f  noise  can  be 
observed  from  a  single  carrier.  A  single  carrier  will  only  provide  a  pulse  in  the 
detector.  Many  carriers  are  needed  to  produce  the  i/f  noise  effect,  just  as  in  the 
case  of  electron  diffraction  patterns.  While  incoming  carriers  may  have  been 
Poisson  distributed,  the  scattered  beam  will  exhibit  super  -  PDissonian  statistics 
or  bunching  due  to  quantum  i/f  noise.  The  quantum  1/f  effect  is  thus  a  two  - 
particle  effect,  best  described  through  the  two  -  particle  wave  function  and  two  - 
particle  correlation  function. 

Let  us  estimate  the  magnitude  of  the  quantum  1/f  effect  by  starting  with  the 
classical  (Larmor)  formula  qv^/Sc3  for  the  power  radiated  by  a  particle  of  charge 
q  and  acceleration  v.  The  acceleration  can  be  approximated  by  a  delta  function 
v(t)  =^v«/‘ (t)  whose  Fourier  transform^i'vis  constant.  The  one  -  sided  spectral 
density  of  the  emitted  Bremsstrahlung  power  2q(Av)t/3ci’  is  therefore  also 
constant.  The  number  2q(^v>‘V3hfc3  of  emitted  photons  per  unit  frequency 
interval  is  obtained  by  dividing  with  the  energy  hf  of  one  photon.  The  probability 
amplitude  of  photon  emission  [2q&v)2/3hfc3]i/2  is  given  by  the  square  root  of 
this  photon  number  spectrum,  including  also  a  phase. factor.  The  beat  term  in  the 
probability  density  I'f'l3  is  linear  both  in  this  Bremsstrahlung  amplitude  and  in 
the  non  -  Bremsstrahlung  amplitude.  Its  spectral  density  will  therefore  be  given 
by  the  product  of  the  squared  probability  amplitude  of  photon  emission  with  the 
squared  non  -  Bremsstrahlung  amplitude  which  is  independent  of  f.  The 
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resulung  specti  al  uensity  of  -fractional  probability  density  fluctuations  is 
obtained  by  dividing  with  I'f'l4  and  is  therefore 

lfr4SI  l2(f)  =  3q2(Av)2/3hf  Nc2  =  fclA/fN  =  r2Sj(f),  (1.1) 

where  =  2e2/hc  =  1/137  is  the  fine  structure  constant  andc^A  =  4q2(d"^ji/3hc^  is 
known  as  the  infrared  exponent  in  quantum  field  theory,  and  is  known  as  the 
quantum  1/f  noise  coefficient,  or  Hooge  constant,  in  electrophysics. 

The  spectral  density  of  current  density  fluctuations  is  obtained  by  multiplying 
the  probability  density  fluctuation  spectrum  with  the  squared  velocity  of  the 
outgoing  particles.  When  we  calculate  the  spectral  density  of  fractional 
fluctuations,  the  velocity  simplifies  and  therefore  Eq.  (1.1)  also  gives  the 
fractional  spectrum  of  current  fluctuations  S  j(f),  as  indicated  above,  ^he  quantum 
1/f  noise  contribution  of  each  carrier  is  independent,  and  therefore  the  quantum 
1/f  noise  from  N  carriers  is  N  times  larger;  however,  the  current  j  will  also  be  N 
times  larger,  and  therefore  in  Eq.  (1.1)  a  factor  N  was  included  in  the  denominator 
for  the  case  in  which  the  cross  section  fluctuation  is  observed  on  N  carriers 
simultaneously. 

The  fundamental  fluctuations  of  cross  sections  and  process  rates  are 
reflected  in  various  kinetic  coefficients,  such  as  the  mobility/^fand  the  diffusion 
constant  D,  the  surface  and  bulk  recombination  speeds  s  and  recombination  times 
Zt  the  rate  of  tunneling  and  the  thermal  diffusivity.  Therefore,  the  spectral 
density  of  fractional  fluctuations  in  all  these  coefficients  is  given  also  by  Eq. 

(1).  This  is  true  in  spite  of  the  fact  that  each  carrier  wi  .l  undergo  many 
consecutive  scattering  processes  in  the  diffusion  process.  The  quantum  i/f  noise 
in  the  mobility  and  in  the  diffusion  coefficient  is  practically  the  same  as  the 
quantum  i/f  noise  in  a  single  representative  scattering  event  which  limits  the 
mobility  or  the  diffusion  coefficient. 
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Due  to  the  rapid  relaxation  of  concentration  fluctuations,  the  quantum  3  /f 
fluctuations  of  scattering  cross  sections  will  only  oe  reflected  by  the 
fluctuations  of  the  mobility  and  the  diffusion  constant  of  the  carriers,  and  not  by 
fluctuations  in  the  concentration  of  carriers. 

For  large  devices  the  concept  of  coherent  state  quantum  1/f  noise  was 
introduced1*'  see  ]_« 0  /^£for  details.  In  this  case  the  Hooge  parameter 
och  may  be  written 


«X  H  =  <*H>coh  =  2 4.4  10' 


where  1/(137)  is  the  fine  structure  constant.  This  is  of  the  same  order  of 
magnitude  as  the  empirical  value  «j_;  =  2  10"^  that  Hooge  found  for  long  devices. 
It  is  therefore  proposed  that  Hooge's  emprical  value  fcr°()_j  is  due  to  coherent 
state  quantum  1/f  noise,  so  that  it  has  a  very  fundamental  origin. 

For  small  devices  <e.g.(  of  sire  L  <  \QUm)  we  apply  conventional,  or 
incoherent^-1*,  quantum  1/f  noise  which  is  just  the  cross  section  fluctuation 
introduced  above  in  Eq.  (1.1).  In  that  case<*H  may  be  written 


o< h  *  ^n’incoh  =  <4w/3F)C(av)<-/(c4-)D, 


where  av  is  the  change  in  the  velocity  of  the  carriers  in  the  interaction  process 
considered.  This  expression  holds  for  any  1/f  noise  source  describable  by 
fluctuating  cross  sections.  Since  usually  (av^/c^X^Tl ,  except  for  carriers  with  a 
very  small  effective  mass,  we  now  have  <  3.3  10“^.  This  may  explain  the  low 
values  of  *ln  the  range  of°(^  =  10-^  -  10-^)  for  very  small  devices.  In 
between  one  can  introduce  a  parameter  s  =  f(L/L0)  where  Lp  is  a  characteristic 
size  and  write 


*H  *  <Vincohcl/<1  +  5,3  +  <Vcoh<s/(1  +  s)3, 


with  s«  i  -for  L/Lc«i  and  s^>l  tor  L/l_D>^i.  According  to  this  rough 
approximation* S  [_r .  =  10y<m  ter  samples  with  a  concentration  c  of  carriers  of 
lO^cm-^  and  vanes  proocrticnal  to  This  describes  the  transition  -from 

Eq.  (1.2)  to  Ec.  ( i .3)  when  one  goes  to  devices  with  smaller  and  smaller  sices,  as 
we  shall  see  below  in  Eq.  (3.21). 

When  we  apply  Eq.  (i.l)  to  a  certain  device  we  -first  need  to  -find  out  which  are 
the  cross  sections  which  iimit  the  current,  and  then  we  have  to  determine  both 
the  velocity  change  o-f  the  scattered  carriers  and  the  number  N  o-f  carriers 
simultaneously  used  to  test  each  of  these  cross  sections.  The  next  section  deals 
with  the  first  problem,  i.e.,  is  dedicated  to  the  study  of  currents  in  MIS  infrared 
detectors. 

II.  CURRENTS  IN  MIS  DETECTOR  STRUCTURES 

MIS  detectors9are  different  from  photovoltaic  detectors,  because  they  do  not 
contain  a  pn  junction  obtained  by  inhomogeneous  doping,  and  use  an  insulated 
field  plate,  or  gate,  placed  on  top  of  a  homogeneously  doped  narrow  -  bandgap 
semiconductor.  The  gate  is  used  to  control  the  surface  potential,  driving  the 
semiconductor  surface  into  deep  inversion.  The  field  of  the  induced  quasi  -  pn 
junction  obtained  under  the  surface  of  the  gate  in  the  homogeneous  semiconductor 
material  is  used  to  separate  the  carriers  generated  by  photo-electrically  induced 
band  to  band  transitions  just  as  in  a  photovoltaic  device.  The  MIS  infrared 
detector  is  therefore  similar  to  a  capacitively  coupled  photovoltaic  detector, 
without  the  inconvenience  of  inhomogeneous  doping  processes. 

MIS  detectors  are  operated  in  the  pulsed  regime  by  applying  the  gate  potential 
which  creates  the  inversion  under  the  surface  for  a  finite  time  only,  and  by 
applying  subsequently  a  potential  which  flattens  the  energy  bands  near  the 
surface  and  releases  the  carriers  which  had  accumulated  from  photo  electric 


effect  and  dark  current  processes  during  the  proceeding  interval.  The  electrical 
signal  obtained  when  the  carriers  are  released,  i.e.,  during  readout,  is 
prcoorcionai  to  the  number  of  carriers  accumulated,  and  therefore  to  the  total 
current  supplying  the  inverted  volume  under  the  surface  with  minority  carriers 
from  the  bulk  and  from  various  thermal  and  photo  electric  processes  in  the 
depletion,  inversion  and  surface  regions.  This  electrical  signal  is  used  in  order 
to  determine  the  flux  of  infrared  radiation.  For  this  determination,  however,  the 
dark  current  contribution  needs  to  be  subtracted  first. 

The  dark  current  is  the  current  supplying  the  potential  well,  defined  by  the 
inversion  region  under  the  surface,  with  minority  carriers  in  the  absence  of  the 
applied  infrared  flux.  Any  low  -  frequency  fluctuation  in  the  dark  current  will  be 
interpreted  as  a  fluctuation  in  the  major  infrared  flux  signal.  Therefore 
fluctuations  of  the  dark  current  at  frequencies  below  the  readout  frequency  will 
limit  the  performance  of  infrared  detectors.  In  the  pulsed  mode  of  operation 
considered  here,  the  dark  current  is  monitored  only  during  the  inverted  phase, 
when  earners  are  accumulated  in  the  potential  well.  Therefore  the  cross  sections 
and  process  rates  which  control  the  intensity  of  the  dark  current  are  not  observed 
continuously  either.  Nevertheless,  the  quantum  1/f  fluctuations  of  these  cross 
sections  and  process  rates  will  be  the  same  as  if  we  would  have  observed  them 
continuously.  Indeed,  the  changes  in  the  incoming  flux  of  electrons  testing  all 
cross  sections  and  process  rates  in  the  semiconductor  is  only  slightly  affected  by 
the  applied  gate  voltages,  and  is  present  also  in  thermal  equilibrium.  This 
independence  of  1/f  noise  on  the  continuous  or  discontinuous  character  of  any 
applied  bias  has  been  experimentally  verified during  the  last  2  decades,  and 
has  been  found  to  be  in  agreement  with  the  interpretation  of  I/f  noise  in  terms 
of  fundamental  resistance  fluctuations.  Although  the  experimental  verification 
was  performed  on  fluctuations  in  conduction  only,  from  the  concept  of  Quantum 
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1/f  noise  we  Know  -that  the  similarity  of  quantum  l/'f  noise  in  the  continuous  and 
pulsed  regimes  should  be  also  t''ue  -for  quantum  i/f  -fluctuations  in  recombination 
cross  sections  and  tunneling  rates. 

The  dark  current  has  to  be  subtracted  -from  the  total  current  in  a  (HgCd)Te  MIS 
device  to  yield  the  photocurrent.  Therefore,  the  minority  -  carrier  dark  current  is 
the  single  most  important  parameter  tor  the  operation  o-f  MIS  devices  as 

g 

detectors  infrared  radiation  .  This  applies  both  to  operation-  of  MIS  devices  in 
the  thermal  equilibrium  mode,  in  which  the  dark  current  determines  the  MIS  diode 
impedance*  and  to  operation  in  the  dynamic,  or  integrating  mode,  in  which  the  gate 
voltage  is  pulsed,  and  in  which  the  minority  -  carrier  dark  current  determines  the 
storage  time  o-f  the  device.  The  main  component  o-f  the  dark  current  in  narrow 
bandgap  HgCdTe  is  the  tunneling  current  via  bandgap  states9,  which  can  also  be 
considered  as  an  electric  breakdown  e-f-fect.  In  general,  the  tunneling  current 
occurs  both  through  band  to  band  transitions  and  through  intermediary  states. 

The  band  to  band  tunneling  current  through  a  simple  triangular  barrier  is 

Jtb  =  (q3E<^/4X3h2)<2m*/Eg)1/2e>:pC4(2m*)Eg3/2/3qhED,  (2.1) 

where  E  is  the  electric  -field  associated  with  the  barrier,  and  Eg  is  the  bandgap. 
The  electric  -field  can  be  approximated  by  the  electric  -field  at  the  semiconductor 
surface 


Es~  (2qn0  <})/££«,> 1/2,  (2.2) 

where  represents  the  empty  well  surface  potential,  and  n0  is  the  doping 
concentration.  Substituting  this  value  into  Eq.  (2),  with  m*/mc,  =  7  10_2EC,  we 
obtain9 

Jtb  =  i0"2no1/2<f>s3/2e>:p[-  4.3  1010Eg2/(no<fe)1/2:  A/cm2,  (2.3) 
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where  n0  is  in  cm-^  and  Eg,1^  in  volts.  Therefore  the  tunneling  current  is 
strongly  dependent  on  the  bandgap  and  also  depends  on  the  doping  concentration 
and  the  surface  potential. 

Experimental  values  of  the  tunneling  current  are  usually  larger  than  Eq.  (2.3) 
because  of  the  additional  effect  of  tunneling  via  bandgap  states.  This  effect  is 
particularly  important  in  n  -  type  devices*®.  Indeed,  in  n  -  type  devices  the 
applied  gate  voltage  is  negative  in  order  to  produce  depletion  at  the  surface.  The 
energy  bands  are  therefore  curved  upwards  at  the  surface,  and  transitions  of 
electrons  from  the  valence  band  to  Shockley  -  Read  (SR)  states  at  the  middle  of 
the  bandgap,  as  well  as  the  subsequent  transitions  from  these  states  to  the 
conduction  band  are  facilitated  by  the  presence  of  many  defects  right  at  the 
surface  of  the  semiconductor.  In  p  -  type  devices  the  similar  indirect  tunneling 
processes  occur  father  away  from  the  surface,  because  in  this  case  the  bands  are 
curved  downwards  at  the  surface,  and  transitions  of  electrons  from  the  valence 
band  to  the  centers  at  the  middle  of  the  bandgap,  as  well  as  the  transitions  from 
the  centers  to  the  conduction  band  well  at  the  surface,  occur  right  where  the 
curvature  begins,  i.e.,  further  away  from  the  surface.  We  conclude  that  in  p  - 
type  devices  there  will  be  fewer  SR  centers  active  in  indirect  tunneling,  and 
therefore  the  tunneling  current  Jc  via  SR  centers  at  a  given  teperature  and  a 
given  applied  gate  voltage  will  be  smaller.  The  tunneling  current  will  be  further 
reduced  in  p  -  type  devices  due  to  the  lower  density  of  states  present  in  the  j 

surface  potential  well  due  to  quantisation  of  the  motion  of  the  electrons  in  the 
potential  well  at  the  surface.  The  reduced  values  of  the  dark  current  in  p  -  type  > 
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devices  correspond  to  higher  values  of  the  breakdown  field  in  these  devices.  The  '  ■ 

best  measured  value*®  of  the  breakdown  field  in  10  cutoff  p  -  type  devices  is  • 
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in  e>:e ss  of  1.0  V/^m,  whereas  that  for  n  -  type  material  of  similar  bulk  defect  ; 
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quality  is  0.5  V !jx m.  On  the  other  hand,  the  minority  carriers  diffusion  current  is  • 

t 

J 


> 


w 

V 

V 


laroer  in  p  -  type  devices  due  to  the  smaller  mass  and  higher  diffusion  constant 
cf  electrons  compared  to  holes.  The  advantage  of  p  -  type  devices  is  therefore 
considerable  only  in  the  case  of  very  narrow  bandaap  and  very  long  cutoff 
wavelengths.  We  shall  therefore  consider  both  the  case  of  p  -  type  and  n  -  type 
devices.  The  large  diffusion  current  present  in  p  -  type  devices  corresponds  to 
the  large  value  of  the  diffusion  length  of  electrons  and  can  be  reduced  by  thinning 
the  devicei  i.e.i  by  reducing  its  thickness  well  below  the  diffusion  length. 

In  general  the  dark  current  Jj  can  be  written  in  the  form 

“d  =  Jdiff  +  Jdeo  +  Js  +  Jtb  +  Jtc  +  Jtsc  +  Jb  +  t2'4> 

The  eight  terms  on  the  right  hand  side  correspond  tQ  minority  carrier  diffusion 
from  the  bulk,  generation  from  SR  centers  in  the  depletion  region,  generation 
from  SR  centers  ai  the  surface,  band  to  band  tunneling,  tunneling  via  SR  centers, 
tunneling  via  surface  centers,  recombination  on  the  back  surface  and  photoelectric 
generation  by  the  thermal  radiation  background  flux<j>p.  We  shall  give  the 
formulae*®  which  determine  these  terms  below,  also  including  an  e>;ample  of  their 
calculation  in  a  p  -  type  device 

Jdif  =  <qni2/p0)CkT^qTn:i1/2  =  (1.6  iO~19C  36  iO24  cm'6/i015cm"3)C0.01  V 
1.5  i05(cm2/Vs)/i0_6sDi/2  =  2  10“4A^cm2,  (2.5) 

Jdepl  =  qniW/2To  =  CU.6  10~19C  6  10l2cm-3  2  10-4cm  /(2  iO“6s): 

*  10-4A/cm2  (2.6) 

Jc  =  Jq  =  qn^s/2  =  C<i .6  iO_i9C  6  1012cm-3  20cm/sD  =  l.S  10_5A/cm2,  (2.7) 

Jtc  =  10"13(Nr4/EQ)  e:cpC-5.3  106Eg2/E3  =  10_i3(4  1012cm_30.2 V/0.062V) 
e>:pC-5.3  i06(0.062V)2/3000 V/cmD  =  8  10-4A/cm2,  (2.8) 


and  Jtb  was  given  by  Eq.  (2.3),  yielding  iQ~^A/cm2.  Here  n,  =  6  10*2cm-2  1S  -the 
intrinsic  carrier  concentration,  and  the  concentration  of  holes  was  taken  to  be  pD 
s  The  mobiiity^n*  i.5  10***cni2/Vs  as  well  as  the  life  time  10-£*s 

or  the  minority  carriers  have  to  be  replaced  by^p  and  Zf>  tor  the  case  of  n  -  type 
devices.  The  surface  recombination  speed  was  denoted  by  s  =  40cm/s,  and  the 
concentration  of  intermediary  states  effective  for  tunneling  was  denoted  by  Nr  * 
n  10*2cm-^.  The  bandcap  considered  was  E_  =  0.065&V,  the  surface  Dotentiaifs  = 
0.2V  and  the  electric  field  below  the  surface  E  =  BOOOV/cm.  All  numerical  values 
have  oeen  included  only  as  an  example  and  are  not  characteristic  of  a  particular 
device.  The  numerical  factors  in  Eq.  (2.6)  correspono  to  p  -  type  HgCdTe  with  a 
20 n  m  cutoff  wavelength  and  were  taken  from  the  paper  of  Kinch  and  Beck^.  For 
n  -  type  devices  we  also  need  to  include  tunneling  via  surface  states.  If  the 
density  of  fast  surface  states  is  denoted  by  N+s  =  !012  cm“2\/-1,  the  current 
generated  from  tunneling  via  a  uniform  density  of  fast  surface  states  across  the 
bandgap  will  be  given  by^ 

Jtsc  =  -qy3lNfskT/2q  =  i.75  1<T4  A/cm2.  (2.9) 

The  sum  of  the  first  seven  currents  on  the  right  hand  side  of  Eq.  (2.4)  must  be 
smaller  than  the  eighth  which  corresponds  to  the  thermal  radiation  background 
flux,  for  background  limited  (BLIP)  operation.  Although  not  all  terms  in  the  dark 
current  are  of  importance,  we  still  retain  them  at  this  point,  because  their 
quantum  1/f  noise  may  be  quite  significant,  even  if  the  corresponding  current  is 
negligible.  We  shall  now  proceed  with  the  calculation  of  quantum  i/f  noise 
contributions  from  all  these  currents. 
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III.  QUANTUM  1/f  NOISE  SOURCES 


1.  1  / -f  Noise  in  the  Diffusion  Current 


The  diffusion  limited  dark  current  will  exhibit  1  /-f  noise  due  to 
conventional  quantum  1/t  -fluctuations  in  the  scattering  cross  sections  o-f  the 
carriers  due  to  phonons  and  impurities.  We  apply  the  -fundamental  -formula  civen 
by  Eq.  <i.i)  -for  an  individual  scattering  process  in  which  the  velocity  change  av  is 
given  by  the  thermal  energy  o-f  the  carriers,  with  the  assumption  that  th? 
collisions  are  perfectly  randomizing  collisions.  If  the  velocity  "v^s  rotated  by  an 
angle  0  in  an  elastic  collision,  the  velocity  change  isl^?)*  2v  sin(0/2).  Averaging 
over  all  anoles  and  velocities,  we  obtain 


6%2  =  4v2  sin2  0/2)  =  2v2. 


and  therefore  from  Eq.  (i.i)  we  get  in  thermal  equilibrium  at  the  temperature  T 
the  i/f  noise  coefficient 


=  <4ot/370<6kT/m*cz), 


where  we  have  assumed  a  Maxwell  distribution  of  velocities.  For  Hgj_xCd>:Te 
with  :<  «=  0.2  we  have  mn*  -  0.008mo  and  for  x  =  0.3  we  have  mn*  =  0.02mo. 
Therefore  we  obtain©^  *  2  in  the  first  case  and^^  =  7.5  10-®  in  the  second 


For  the  case  of  Umklapp  scattering,  which  occurs  in  semiconductors  only  to  a 
limited  extent  due  to  the  relatively  small  number  of  high  momentum  phonons 
available  at  the  temperature  T,  the  momentum  change  of  the  electron  is  given  by 
the  smallest  reciprocal  lattice  vector,  and  therefore^  v  =  h/am*.  We  therefore 
obtain  the  quantum  1/f  noise  coefficient 


=  <4<*/37D(h/m*ac)2, 
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which  is  much  larger  than  Eq.  <3.2)«  but  has  to  be  multiplied  with  a  negative 
exponential  which  describes  the  scarcity  of  phonons  with  momentum  of  the  order 
of  a  reciprocal  lattice  vector.  CThe  negative  exponential  could  be  included 

in  the  current  weight  factor  which  will  be  defined  below  in  Eq.  (4$,  but  we.preter 
to  include  it  here  already].  Combining  Eqs.  (3.2)  and  (3.3),  we  obtain  -for 
conventional  1  /f  noise  in  the  mobility  and  diffusion  coefficients 

<*H  =  (4«/3^r)C(6kT/m#c2)  +  (h/m*ac)2exp(-^/T)3,  (3.4) 

where  &  is  about  half  the  Debye  temperature  for  simple  metals,  but  may  be 
higher,  of  the  order  of  the  Debye  temperature,  for  semiconductors. 

The  quantum  l/'f  noise  considered  so  far  is  known  as  conventional  quantum  1/f 
noise,  and  affects  cross  sections  and  process  rates.  In  sufficiently  large 
semiconductors  samples  we  expect  a  larger  form  of  quantum  i/f  noise,  described 
in  Appendix  A  and  known  as  coherent  state  quantum  1/f  noise.  For  this  type  the 
1/f  noise  coefficient  is  given  by 

«coh  =  2oc/jr  =  4.6  10~3.  (3.5) 

The  values  of  the  quantum  1/f  noise  coefficient  given  by  Eqs.  (3.1)  -  (3.5)  can 
be  used  to  calculate  the  quantum  i/f  noise  which  affects  the  various  currents 
listed  in  Eq.  (2.4).  We  first  consider  the  case  of  the  dark  diffusion  current  of 
electrons  from  the  bulk  through  the  surface  barrier  in  a  p  -type  MIS  device, 
similar  to  diffusion  in  a  n+p  junction,  because  in  both  cases  the  current  is 
determined  by  the  diffusion  of  electrons  which  are  minority  carriers,  against  the 
built  -  in  field  of  a  Boltzmann  potential  barrier  into  the  surface  well,  and  by  the 
thermal  generation  of  carriers  there.  We  start  with  the  derivation  of  the 
mobility  fluctuation  part  of  quantum  i/f  noise  in  a  n+-p  diode.  For  the  MIS 
barrier,  just  as  for  a  diffusion  limitted  n+-p  junction,  the  current  is  controlled  by 


-6.0- 


dilution  of  electrons  into  the  p  -  region  over  a  distance  of  the  order  of  the 
diffusion  length  L  =  (Dn'T(_() 1  which  is  usually  shorter  than  the  length  wp  of  the 
p  -  region.  If  Nix)  is  the  number  of  electrons  per  unit  length  and  Dp  their 
diffusion  constantf  the  electron  current  at  x  is 

Ind  =  -  eDndN/dx,  (3.6) 

where  we  have  assumed  a  planar  junction  and  taken  the  origin  x  =  0  in  the  junction 
plane.  Diffusion  constant  fluctuations)  given  by  kT/e  times  the  mobility 
fluctuations*  will  lead  to  local  current  fluctuations  in  the  interval^  x 

Jlnd(x,t)  =  I_d4nri(x,t>/Dn.  (3.7) 

The  normalized  weight  with  which  these  local  fluctuations  representative  of  the 
interval^  contribute  to  the  total  current  Id  through  the  diode  at  x  =  0  is 
determined  by  the  appropriate  Green  function  and  can  be  shown  tD  be 
(Ax/L)exp(-x/L)  for  wp/L  >>  i.  Therefore  the  contribution  of  the  sectioruPx  is 

ckJd<>:,t)  =  (Ax/L)exp<-x/L)IndcfDn(>:,t)/Dn,  (3.8) 

with  the  spectral  density 

S^idbef)  =  (4x/L)2  exp(-2x/L)  Ind2  SDn(x*f>/Dn2.  (3.9) 

For  mobility  and  diffusion  fluctuations  the  fractional  spectral  density  is  given  by 
c^Hnd^M^’  wher^-Hnd  *s  determined  from  quantum  i/f  theory  according  to  Eos. 
(3.1)  -(3.5).  With  Eq.  (3.6)  we  obtain  then 

S^id'Xif)  =  &x/L2)  exp(-2x/L)  (eDpdM/dx)^  pnd/fN.  (3.10) 

The  electrons  are  distributed  according  to  the  solution  of  the  diffusion  equation, 

i.e. 


v.o.-.'vy 


N(x>  >=  C N (0)  -NpD  e>:p(-x/L)  +Np;  dN/dx  =  -CCN(O)  -Kp3/L>  exp(-x/L>.  (3.11) 

Substituting  irrto  Eq.  (3.10)  and  simply  summing  over  the  uncorrelateo 
contributions  of  all  intervals  Ax,  we  obtain 

A 

sIdw  =c<Hnd(eDn/L2)2|  CN<0)  “  Np32e“4>:/Ld>:  /<CN(0)  -  K_3e">!/L  +ND>  (3.12) 

We  note  that  (eDn/Lt)2  =  <e/Tn)2.  With  the  expression  of  the  saturation  current 
Iq  -  e(Dn/,cn)1  ^‘•Np  and  of  the  current  I  =  IgCexpleV/kT)  -  13  we  can  carry  out  the 
integration 

I 

sId(f)  s  wHnd(eI/'fTn,/a2u3du  /(*u  +  1} 

s  aHnd<el/fV  £F(a) -F(a  e“w /L)3 -*Hnd(eI/f^n)aw/C(a  +  1)L3,  (3.J3) 

the  last  form  being  for  wp  <<  L.  Here  we  have  introduced  the  notations 
u  =  exp(-x/L),  a  =  exp(eV/KT)  -  1, 

F(a)  =  1/3  -  l/2a  +  1/a2  -  (l/a3)ln(i+a>.  (3.14) 

For  wp  >  >  L  we  have  F(o>  =  0,  and  the  second  term  in  rectangular  brackets  drops 
out  in  Eq.  (3.13). 

Eq.  (3.13)  only  contains  the  fluctuations  in  the  mobility  and  the  diffusion 
constant.  In  a  similar  way  we  calculate  the  quantum  1/f  fluctuations  of  the 
recombination  rate  in  the  bulk  of  the  p  -  region.  We  have  for  the  recombination 
current  Alp(x)  in  a  section  Ax,  if  N'(:<)  is  the  excess  carrier  density, 

AIR(x)  =  eN'(x)Ax/Tn  .  (3.15) 

Putting  Cn  =  1/Tn  and  bearing  in  mind  thatTn,  and  hence  Cn,  fluctuates,  we  have 


for  the  section  Ax, 
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Scr.'-'+'-'W,-  =^hpr/'fAN 
■  ■  that,  since  Mo;)  = 


<3.U) 


3^1 -o: ,  4)  =  AI-^oOE^- j^pp/'SAND 

*ot  ^ppCe^CN'oMZ^  /tin“N(>;)}Ax 


(3.17) 


whe^e  !!'(>:)  =  CMiO'  -  Np3e>:p(->: / L)  end  No;)  =  i! ' (>: )  +  Np  as  before. 

It  is  easily  snown  that  the  -fluctuating  current^ !(>;,  t)  at  the  junction  is 


Jills,  t)  =Ja  Ip(>: i  t)e:<p<-:>- /L) 


(3.13) 


f 

Sj (t;  =0(Hnp\e2NpL/-frn2)  )  {CN'(x)32/NDM(>;)3e>;p<-2>:/L)d(x/L)  = 

Hnr;eID.-fTni  f  [a2u3/(au  +  l)Ddu  =Bt,^nrC9l/-frr(jLF(a)  -  F(ae_W//L)3  (3.19) 

ep^r/L 

where  F(a),  a,  and  u  have  the  same  meaning  as  before. 

We  can  use  the  similarity  o-f  the  quantum  i/f  noise  results  for  diffusion 
current  fluctuations  caused  by  mobility  fluctuations  and  by  recombination  speed 
fluctuations  in  order  to  combine  both  into  a  single  formula 

Sj(f)  =(d  Hd  +^Hnr)CeI/fi:n3CF(a)  -  F(ae~w/L)3.  (3.20) 

In  the  limit  of  very  short  devices  <wp  «  L)  the  last  factor  becomes  aw/C(a  +  1)Lj, 
and  in  the  limit  of  long  MIS  devices  <wp  >>  L)  it  simply  becomes  F(a).  In  addition 
we  have  a  current  noise  contribution  sJb  from  the  quantum  1/f  fluctuation  of  the 
recombination  speed  s  on  the  back  surface. 

Sc  far  we  have  considered  only  conventional  Quantum  1/f  noise  which  is 
applicable  to  sufficiently  small  devices.  In  general,  however,  we  must 
interpolate  between  conventional  and  coherent  quantum  1/f  noise,  according  to 
the  relation 


=  Cl/<1  *  s)3C2  A/f3  +  Cs/d  +  s)2Z2oi/7ir j, 


—  6  3“ 


where  s  =  Em  /E^  =  2e“N'/m*c“,  and  N'  =  nS  is  the  number  o f  the  carriers  per  unit 
length  of  the  device,  s  represents  the  ratio  between  the  magnetic  energy  per  unit 
length  and  the  kinetic  energy  per  unit  length  ct  the  device.  The  quantity  e"-/m*c- 
*  ^rn/m*  can  be  calculated  in  terms  of  the  classical  radius  of  the  electrpn  rc  = 
2.84  10~13  cm.  Then  we  obtain  s  =  2  r.-N'm/m*,  i.e.,  the  parameter  s  represents 
twice  the  number  ot  carriers  present  in  a  length  o-f  the  device  equal  to  the 
classical  radius  ot  the  electron.  We  must  compare  s  with/TA,  and  it  s  <<7t  a  we 
apoly  conventional  quantum  i/'-f  noise ,  whereas  -for  s  >>A  A  we  have  to  apply 
coherent  quantum  i/f  noise.  In  general  the  approximate  -formula  o-f  Ec.  (3.21) 
must  be  used  -for  the  transition  region. 

The  dimensionless  parameter  s  is  easy  to  calculate  in  any  practical  case.  For 
instance  in  the  case  o-f  a  MIS  device  o-f  area  50yU.tr>  x  50/tm  with  a  concentration  o-f 
carriers  o-f  10^  crrT^  we  obtain  N*  =  2.5  10^ /cm,  and  with  m/m*  =  50  we  obtain  s 
=  7.  On  the  other  hand,  we  can  estimate  A  tor  conventional  quantum  1/f  noise  and 
we  will  certainly  find  A  «  i,  because  the  velocity  change  of  the  carriers  must  be 
much  smaller  than  the  speed  of  light.  Therefore,  in  this  case  we  must  apply 
coherent  quantum  i/f  noise,  because  s  »  7TA.  Consequently,  in  Eq.  (3.20)  we  must 


^  Hd  +0^  Hnr  =o6coh  =  A-6  10_3 


(3.22) 


The  coherent  state  quantum  i/f  noise  coefficient  thus  replaces  the  total 
conventional  Hooge  parameter. 


2.  1/f  Noise  of  the  Recombination  Current  Generated  in  the  Depletion  Region 


The  quantum  1/f  noise  of  the  recombination  current  thermally  generated  in  the 
depletion  region  arises  from  quantum  1/f  fluctuations  of  the  bulk  recombination 


s.’r.sx.t.s 


rates  in  the  depletion  region.  The  difference  between  the  recombination  rate  R 
and  the  generation  G  is  given  by 


m 


=  -  G  =  Con  -  n-jD/Cin  +  n^r  +  <p  «■  Pj)rno] 


(3.23) 


v/here  nj  and  pj  are  electron  and  hole  densities  when  the  Fermi  level  lies,  at  the 


trap  level.  It  the  trap  level  lies  at  the  intrinsic  level,  nj  =  pj  *  n^.  Moreover. 


cpo  andfno  are  time  constants  tor  electrons  and  holes.  If  A  is  the  cross  - 


sectional  area  of  the  junction,  the  current  is 


<u*r  vs 

1  =  J e(R  -  G)A  d;:  =  e  j  Cpn  -  n^jl/nn  +  n,)rpD  +  <p  +  n^)tnoDA  d>:  (3.24) 

0  J 


where  w  is  the  width  of  the  space  -  charge  region  and  the  trap  level  is  assumed 


to  lie  at  the  intrinstic  level. 


We  now  turn  to  the  a  -  r  noise.  The  time  constants Trpo  and  T  nD  fluctuate  in  a 


1/f  fashion  and  this  produces  the  quantum  i/f  contribution  to  g  -  r  noise.  We  now 


write 


^”no  =  l^n 


po  =  ^Cp 


(3.25) 


where  Cn  and  Cp  are  the  generation  (or  combination)  rates  for  a  single  electron 


and  for  a  single  hole,  respectively.  Consequently 


^  L  no/^no  =  ~^cn/cn) 


tTpo/Zpo=  -<Jcp/Cp) 


(3.26) 


We  now  aoply  this  knowledge  to  Eq.  (3.23)  and  observe  that 


<5(R  -  G)  =  CR(:<)  -  GO:)] 


U<n  +ni)rM>tfCD/CJ3  +C(p  +Pl>rnD/cn/Cn)3>/C(n  +n jT  +  -ip  +ni)Z’nc3.  (3.27) 


m 


km; 

m 

g 

LrpiLi 

1 


I 


i 
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S  i  =  e  I  J'CR(X)  -  G(x)  JAd::> 


(3 .25) 


the  noise  is 


'K  * 

=  e2/ j  lR<::>  -  G(x)3AC(R(x')  -  G(x')3A 


CC(n  +  nt)2t  pC,25cp<«*  «'*  *)/Cp2  +  <P  +  ni,2rnoisCn<>:’  :</»  f)/Cn"3 


/Cln  +  n^T  +  ip  + 


(3.29) 


since  5  C_  and  6  are  independent. 

u1  »  I 


We  now  observe  that 


pn  -  n2^  =  n^jCe>:p(e  V/2KT)  -  1  3te>:p(e  V/2KT)  +  13 


(3.30) 


and  that  the  integrand  in  Eq.  (3.29)  has  an  appreciable  value  only  i-f  p~  n  -  ni 


exp»e V/2KT).  By  substituting  n  +  ^  =  p  +  r.j  =  niCe>:p<eV/2k’0  +  13  we  de-fine  an 


e-f-fective  width  such  that 


eA  J  C(pn  -  n^dxD/Ctn  +  n^r^  +  <P  +  nii2Yio3 
=  eA((pn  -  n2j)/n^Cexp(e V/2KT)  +  1  3Hwe.ff/(E'po  no)3  (3.31) 


We  may  thus  write 


I  *  lgr  =  eAwe^niCexp(eV/2KT)  -  13/(Zpo  +2Tno) 
=  CeN,„/£po+rnomanh  eV/2KT 


(3.32) 


where  =  AweHnjCexp(e V/2KT)  +  1 3  is  the  e-f-fective  number  o-f  hole  -  electron 
pairs  taking  part  in  the  conduction  and  noise  processes.  This  equation  is  exact 


but  not  very  use-ful  since  it  contains  the  unknown  parameter  w^. 


We  now  turn  to  Eq.  (3.29)  and  observe  that 


/  ,V 


and 


SCpls,  x',  -f)/C:D  =  (W  Hp/f)/CR(>:')  +  G(x,:']<rpo  n0)Al/(x'  -  >:)  0.33) 


Hp 


SCn<>:,  x',  f)/C2n  =  {faCHn/f)/CR(xO  *  G(x')3(Tpo  +Tno»A>J(x'  -  >:)  (3.34) 

The  -factor  (rn0  *T  po>  enters  in  because  SCp(>:,  s',  ■f)/C‘Lp  and  SCn<>:,  x',  fl/C-,. 
must  be  independent  o-f  Z p0  and~no  i-f  prn-nj  exp(e V/2KT).  This  yieldsi  i-f  we 
integrate  over  the  $  function 

Sjtf)  =  e2  f  { CR(x)  -  G(:<)32A/(rpo  +  r  no,ER(;:)  *  G(>:);1> 

(C(n  +  nj)2rpo  oi.  Hp/f  +  (p  +  ry^ncV  Kn/f 3/C<n  +  n^r^  +  <p  +  ni>rno32}dx  0.35) 

We  now  observe  that  the  second  factor  in  So.  0.35)  is  practically  a  constant 
as  long  as  p  and  n  are  comparable.  We  may  thus  bring  that  factor  outside  the 
integral  sign  and  write 

C(n  +  ni)1Tpo<fc£Hp^  +  *P  +  ’V^no^Hn^2 
/C(n  +  njltpo  +  <p  +  nplno^2  =d-  H/f  ;  (3.36) 

whereof  ^  is  given  by 

d'  H  =  ^po^po  +  L  no^  Hp  +  L  po  +  ^  no^Hn 


We  thus  have 


Sj <-*>  -I pf  He/f<tno  +~po);i  *  J  ER!>:)  _  G(x!D2A/[R(x>  +  Gtxfld 
"O^He^r'^no4?  po>  2tanh[e  W2KT] 


x 

(3.33) 


We  can  now  prove  Eq.O.38)  in  a  Hooge  -  type  formulation  [ID. 
Here  we  put 


Sj(f)/I2cr  =ot  H/fNeff 


(3.3?) 


But  according  to  Eq.(3.32).  I 


=  <eNe^.f/£/tanh  eV/2KT»  so  that 
Sj(-f)  =oi  HCeIgr/rntar,h  eV/2kT,  X -  <rnD  +Z  po>  (3.40) 

in  agreement  with  Eq.(3.33). 

We  can  also  prove  Eq.(3.3S)  -from  the  -following  consideration.  We  write  I_p  = 
(No^f^D  tanh  eV/2kT,  where  M  =  e /^fluctuates.  We  then  have 

SAl<-f)/d  I2  =d  H/-f  or  SAl(-f)  =CC  H/-f(e/2D/4l  (3.41) 

so  that)  since  the  hole  -  electron  pairs  are  independent 

S  :(f)  =  Ng^S^jW  =  (eNeHdI/T'Z;^CH  =  C(eIgr/fT6CH3tanh  ev/2kT  (3.42) 

whereof  ^  is  given  by  Eq.<3.37)  and  T  s.(Tno  +Z  pC,). 

The  last  two  approaches  are  easily  extended  to  other  cases;  the  method  works 
as  long  as  a  time  constant  X  and  an  can  be  de-fined. 

We  -finally  evaluated.  andd  -from  quantum  i/f  noise  considerations  C23» 

m. 

d  Hn  3  <4d/3X)<Av2n/c2)  s  4o!/3;rC[2ea(Vdif  -  V)  +  3kT3/m*nc2>  (3.43) 

ot  Hp  =  <4d/370G4v2p/c2)  =  4d/37TlC2e(l  -  a)(Vdif  -  V)  +  3kT3/m*pc2l  (3.44) 

and  as  a  consequence  (see  Eq.(3.37)  and 

^  H  =  <4d/370(C2e<Vdi<  -  V)  +  6kT3/C(m*n>1/2  +  <m*p)1/232c2  (3.45) 

The  problem  has  hereby  been  solved.  Note  that  in  Hgj_>;CD>;Te  with  ::  =  0.3.  m*n 
=  0.02  m,  m*p  =  0.55  m,  so  that  C(m*n>1/2  +  <m*p)1/232  =  0.76?  m,  very  much  larger 


3.  Noise  in  "the  Surface  Recombination  Current 

The  surface  recombination  currnet  is  a  dark  current  component  originating 
from  surface  states  at  the  interface  between  the  surface  passivation  layer  and 
the  bulk.  The  recombination  cross  sections  of  these  states  will  exhibit  quantum 
i/f  noise.  The  quantum  i/f  noise  coefficient^  ^  reflects  the  velocity  change  of 
carriers  involved  in  this  process 

oi  H  =  <4oi/37n<2/m*c2)C3kT/2  +  eU/2  +  Q.leVl,  (3.46) 

where  U  is  the  surface  potential  jump  present  between  the  surf  ace. passivation 
layer  and  the  bulk  even  if  no  gate  voltage  is  applied,  and  V  is  the  applied  gate 
voltage  which  we  take  with  a  coefficient  less  than  unity,  here  for  example  with  a 
weight  of  0.1 . 

The  calculation  of  surface  recombination  quantum  i/f  noise  is  similar  to  the 
calculation  of  quantum  i/f  noise  from  recombination  in  the  space  charge  region. 
However,  in  this  case  the  cross  sections  are  not  distributed  over  the  width  of  the 
Junction,  but  rather  are  concentrated  at  the  surface  which  is  caractenced  by  the 
surface  potential  Therefore  we  can  write  an  expression  of  the  form 

SJs(f)  HeJstanh  x  /f<Tno  +T  po>e>;p(e  <J>/kT).  (3.47) 

We  note  that  the  fabrication  process  of  MIS  structure  introduces  less  bulk 
defects  than  the  fabrication  of  photovoltaic  devices.  Nevertheless,  the 
fabrication  of  MIS  devices  introduces  some  defects  in  the  bulk  layer  located  right 
under  the  surface.  These  defects  will  manifest  themselves  through  a  contribution 


to  indirect  tunneling. 


4.  Quantum  l/f  Noise  in  the  Tunneling  Rate 

In  the  case  of  tunneling  -from  a  surface  accumulation  layer  to  the  bulk,  the 
velocity  change  o-f  the  carriers  will  lead  us  -from  the  thermal  velocity  o-f  the 
carrier  on  one  side  o-f  the  barrier  to  the  thermal  velocity  o-f  a  carrier  on  the.  other 
side  o-f  the  barrier,  i-f  band  to  band  tunneling  is  considered.  I-f,  however,  the 
tunneling  goes  via  intermediate  states  located  in  the  banrigap,  the  velocity  is 
zero  as  long  as  the  carrier  is  stationary  in  the  intermediate  state.  We  can 
therefore  write  the  i/f  noise  coefficient 

Of =  (4a/3'T)6kT/m#c2  (3.48) 

for  band  to  band  tunneling,  and 

*  <4*/37r)3k7/m*c2  (3.4?) 

for  tunneling  via  intermediate  states  in  the  bandgap,  where  we  have  considered 
the  average  squared  velocity  change  two  times  smaller.  The  effective  mass  is  the 
mass  of  the  minority  carriers  in  the  bulk  material. 

For  band  to  band  tunneling  from  a  surface  inversion  layer  to  the  bulk,  the 
velocity  change  of  carriers  corresponds  to  an  energy  difference  of  the  order  of 
the  bandgap  Eg  plus  an  energy  differnce  of  the  order  of  the  thermal  energy  3kT/2, 
provided  we  are  dealing  with  deep  inversion,  as  used  in  practical  MIS  devices  in 
the  pulsed  mode  of  operation.  This  yields  the  quantum  i/f  coefficient 

CK  -  (4of/37T)(Eg  +  3kT/2)/m*c2  (3.50) 

for  sand  to  band  tunneling.  For  tunneling  via  intermediate  states  in  the  bancgap 
the  corresponding  energy  differnce  will  be  smaller,  and  therefore  we  replace  EQ 


-  ~70  - 


\ 

• 

0^=  <4«i/3‘ir)(Eg/2  +  3KT/2)m*c2  (3.51)  | 

f 

I 

This  relation  is  applicable  both  it  the  intermediate  states  are  located  in  the  i 

t 

depletion  region  or  at  the  surface.  ' 

i 

In  the  last  four  equations  we  did  not  divide  by  the  number  of  carriers 

i 

simultaneously  involved  in  the  tunneling  process,  because  this  number  is  less 
than  unity  for  practical  tunneling  currents.  Whenever  a  cross  section  or  a 
process  rate  is  tested  with  one  electron  or  less  than  one  electron  at  a  time,  the 
effective  number  of  electrons  in  the  denominator  of  the  quantum  1/f  formula  must 
be  replaced  by  unity.  Let  us  calculate  the  average  number  of  carriers 
simultaneously  present  in  the  tunneling  process  at  any  time.  The  tunneling 
process  occurs  over  a  distance  d  =  En/eE,  and  the  speed  v  of  the  carriers  will  be 
of  the  order  of  the  thermal  speed  in  the  case  Df  an  accumulation  layer,  and  of  the 
order  of  the  bandgap  energy  in  the  case  of  a  surface  inversion  layer.  Dividing  d 
by  v,  we  obtain  a  tunneling  time  t  *  6'iO“12s  for  accumulation  layers  and  t  =  3- 
10-*2s  for  inversion  layers.  From  Eq.  (2.S)  we  Know  that  the  tunneling  current  is 
of  the  order  of  10-2A/cm2.  Multiplying  this  by  t,  we  obtain  3  -  6‘10-lc>  C,  i.e., 

2000  -  3000  electrons/cm'-  tunneling  simultaneously  in  a  device  of  1  cm'-  area.  In 
a  device  of  dimensions  50  jm  >:  SO^m  =  2.5  10“^  cm2  the  average  number  of 
carriers  in  the  process  of  tunneling  at  any  time  is  therefore  0.05  -  0.075,  and  this 
is  indeed  much  less  than  unity.  Nevertheless,  if  the  area  of  the  device  exceeds  5 
10~*  cm2,  Eqs.  (3.48)  -  (3.51)  require  an  additional  factor  e/tJ*A  which  makes  the 
noise  spectral  density  proportional  to  and  A,  rather  than  to  the  square  of 
these  quantities. 

The  photoelectric  current  will  reflect  the  fluctuations  in  the  number  of 
photons  arriving  from  the  radiation  background.  The  quantum  efficiency  will  not 
exhibit  considerable  quantum  1/f  noise,  because  the  generated  carriers  will  be 
corrected  with  certainty.  Therefore  the  collection  of  photoelectrically  generated 


.i  tf.a  <,i. 


earners  is  not  controled  by  any  cross  section  or  process  rate  affected  by 
considerable  quantum  i/f  noise.  If  we  neglect  the  i/f  noise  generated  in  the 
senes  resistance  of  the  diode*  there  should  be  no  photoelectncally  generated  1/f 
noise  -from  a  short  -  circuited  diode. 

Since  the  various  dark  current  fluctuations  with  1/f  spectrum  are  statistically 
independent,  the  total  1/f  noise  is  simply  obtained  by  summing  all  contributions. 


IV.  1/f  NOISE  LIMITED  PERFORMANCE  OF  MIS  DIODES 


From  Ea.  (2.4)  we  write  the  total  dark  current  fluctuation  in  the  form 

</Jd  =  +cfJcjer,  +4  +/ Jtb  +  cf  Jtc  Jtsc  +  (4A) 

and  the  spectral  density  of  current  fluctuations  will  be  neglecting  4 

3Jd  =  5Jdif  +  sJdep  +  SJs  +  SJtb  +  5Jtc  +  SJtsc-  (4>2) 

Here  we  have  lumped  the  recombination  current  on  the  back  surface  Jb  together 

with  the  surface  recombination  (generation)  current  Jg.  If  we  denote  all  the 
corresponding  spectral  densities  of  fractional  fluctuations  by  a  prime*  S'jj  = 
Sji/J>,  we  obtain 

s'  Jd  =  ('-Tdif/Jd)<'Sdif  +  wdep/,'Jd)*'S,Jdep  +  <Js/‘Jd,i’S,Js 
+  <Jtb/'-Td)2s'jtb  +  <Jtc/Jd)2s'jtc  +  <Jtsc^d)2s'jtsc-  (4'3) 

This  equation  was  obtained  by  dividing  the  previous  equation  through  and 
shows  that  the  biggest  contribution  will  not  necessarily  come  from  the  process 
with  the  highest  fractional  quantum  1/f  noise*  i.  e.»  with  the  highest  1/f  noise 
coefficient.  The  weight  of  each  type  of  noise  is  determined  by  the  corresoonding 
squared  current  ratio. 
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Tne  detectivity  of  infrared  detectors  is  limited  in  general  by  three  types  ot 
noise:  u)  curent  noise  m  the  Getector,  (Id  noise  due  to  background  photons  lohcton 
noise),  nil)  noise  m  the  electronic  system  following  the  detector.  We  shall 
neglect  here  the  background  photon  noise  and  the  noise  in  the  electronic  system. 
The  detectivity  is  defined  as 


B*\X,  t)  =  <AAfli/2/NEP(cmHz1/2/W) 


where  A  is  the  area  of  the  detector)  NEr  the  noise  equivalent  power  denned  as 
the  r .m.s.  optical  signal  o f  wavelenghA  required  to  produce  r.m.s.  noise  voltage 
(current)  equal  to  the  r.m.s.  noise  voltage  (current)  in  a  bandwidth  A  f,  and  f  is  the 
frequency  of  modulation.  The  noise  equivalent  power  NEP  is  aiven  by 

NEP  =  (h V  />£  q)(  5Id(-f)/-f)i/2.  (4.5) 

Therefore  we  obtain  -for  the  detectivity 

D#(7V,  4)  =  <^q^/hc)  CA/SId(f)J1/2  =  t^qa/hc)  CSJd(t)D*1/2  (4.6) 

We  notice  that  B*(A»  f)  is  proportional  to  A  up  to  the  peak  wavelength  Ag.  For  A  >\c 
we  have”^  =  0  and  thus  D#(t\)  f)  =  0.  By  substituting  our  result  for  Sjd,  we 
obtain  the  general  expression  of  the  detectivity  as  a  -function  of  various 
parameters  o-f  the  MIS  device. 

Let  us  now  evaluate  the  magnitude  of  the  various  dark  current  noise 
contributions.  With  mh*  =  0.55  mD,  mn*  =  0.02  mD,rn  =  10_t55,  Eg  =  0.1  eV,  3kT/2 
-  0.01  eV,  Ns+-  =  1014/Vcm'-)  we  obtain  -for  a  p  -  type  device  w<+n 

S'jdiff  =  *Hnd  +<^  Hnr^e/,^n^diff  ^coh(e1/2/Cf(l(T/n)1/2Np];F(i)/a 
=  (4.6  10_3/4fNp)  4-10“10C1/2/[10‘6s  1 .5*  105(cm2/ Vs)  4-10-21  J31/2 

=  1  .£:'  10  ^cm“/f,  Cor  =  10“^  cm“-/f)  with  0(  =  2’10-^  for  incoherent  noiseD  (4.7) 

Up 


s'jdep  =  C^Hee/^rno  +  po)Jdep3t*nh  ^V/ZkT 
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=  Ci!(j^ee/teAwn<e V/'2KT)3e V/2kT  s°^^p/fAwn1  =  4.6  10-^cm*;/t  (4.S) 

S'j£  =  <4<*/37r>(2/m*c2>i:3kT/2  +  eU/2  +  O.iVe]  Ce  tanhx  /ffTno  +rp0)Jg] 
■  <4«L/37r0.02)(2/500000)C0.025  +  0.5  +  0.5]  le  tanhx  /•feAwniiee V/2kT  -  i)] 

=  7  10~5cm2/f  (4.9) 


S-jifc  =  <4<*/37r>(Eg  +  3kT/2)/m*c2  =  (4/9.5  137  0.02H0.1  i/500000) 
=  3.3  10"acm2/f  (4.10) 

S'j-tc  =  <4«*/3Jr>®g  +  3KT)/2m*c2  =  <4/9.5  137  0.02X0. 12/106) 
a  1.3  10“:::cm2/f  =  S'Jisc  (4.11) 


s,Jdi+-f  was  calculated  In  the  small  bias  limit  -for  wp>>L»  but  wp=  0.25  L  gives  the 
same  result;  the  incoherent  case  with  a  lattice  constant  D-f  0.65  nm  and  $=  320  K 
was  also  listed  above  (because  a  lOyicm  thick  device  is  very  short,  so  it  may  be 
applicable),  and  would  give  1.3  10_10cm2/f  (or  a  n  *  type  device.  Eqs.  (4.8)  - 

-  (ft 

(4.11)  would  be  reduced  only  m*/mrt  =  27.5  times  for  n  -  type  devices.  We  mention 
that  S'j5  has  been  calculated  with  the  inclusion  of  a  term  of  10%  of  the  applied 
gate  voltage  V  into  the  kinetic  energy  of  the  carriers  at  the  surface,  and  that  for 
the  back  surface  recombinaiton  current  this  term  has  to  be  dropped  in  the  similar 
expression  of  S'jb-  However,  we  have  neglected  this  here,  because  the  surface 
recombination  terms  will  not  turn  out  to  be  important,  as  we  will  see  below.  The 
applied  gate  voltage  was  taken  to  be  V  =  5  V.  Using  Eq.  (4.3)  and  the  current 
densities  evaluated  in  Eqs.  (4.10)  -  (4.14)  to  calculate  the  fraction  of  each 
current,  we  obtain 


1  cm-2  f  S'j(f)  =  (20/132)2  1.3  10-6  +  (10/132)2  4.6  10-9 
(0.01/132)2  3.3  10~&  +  (30/132)2  1.3  10“?  +  (17.5/132)2  1.3 
=  3.67  10“®  +  2.6  10"11  +  5.2  10-il  f  1.9  10-16  +  6.61  10“9 


+  (3.6/132)2  7  10“8 

io-s 

+  3.17  10'i0 


+ 


=  4.37  10~s,  or  for  incoherent  1/f  noise,  7.1  10-9  (p)  and  3  10“^ 


(n).  (4.12) 


This  value  can  be  used  in  order  to  estimate  the  detectivity  of  the  device  in  our 


example.  Substituting  into  Eq.  (4.6),  we  obtain  with  a  quantum  efficiency  ^  *  0.7 
and  wavelength  of  ICy^m 


Dfra,  f)  =  qfi/ho  CSJd(f)J"1/:=  C0.7'1.6'10'19C  10"5 
lCi*m/s)J  Cf/(4.37-10-3cm2  i.74,10“*A2/cm4>:)i/2  =  2*107 
for  incoherent  1/f  noise,  5‘107  (p),  and  2.5  10°  <n>.  (4.13) 


m/(6.6-10"34j5  3* 
(cm  Hz^2/w)  f1'1’*-, 


In  conclusion  we  note  that  for  the  relatively  large  devices  which  we  have 
considered,  most  of  the  quantum  1/f  noise  comes  from  fluctuations  m  diffusion 
and  in  the  tunneling  rate  via  impurity  centers  in  the  bandgap.  The  effective  mass 
of  the  earners  is  present  in  the  denominator  of  all  quantum  1/f  noise 
contributions  except  the  coherent  quantum  i/f  fluctuation  present  in  the 
diffusion  current  of  large  devices.  In  smaller  devices  the  diffusion  current  will 
also  be  given  by  the  conventional  quantum  1/f  formula  which  contains  the 
effective  mass  of  the  carriers  in  the  denominator.  For  umKlapp  scattering  the 
mass  of  the  carriers  in  the  denominator  is  even  souared.  Consequently  we  expect 
lower  quantum  l/f  noise  from  n  -  type  devices,  in  which  the  minority  carriers  are 
holes,  particularly  if  the  devices  are  very  small,  e.g.,  below  10/c  m. 


V.  DISCUSSION 

The  transition  from  coherent  state  quantum  i/f  noise  to  conventional  quantum 
1/f  noise  is  particularly  interesting,  and  should  be  studied  experimentally.  This 
is  possible  with  a  sequence  of  devices  of  smaller  and  smaller  size,  and  will  show 
a  considerable  change  in  noise  at  a  size  of  the  order  of  lO^m1-.  The  theory  of 
the  transition  is  not  yet  well  developed.  Therefore,  this  experiment  has 
particular  importance-  we  do  not  know  if  the  parameter  s  is  sufficient  to 


m 
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characterise  the  transission,  and  it  the  parameter  s  should  not  be  replaced  by  a 
power  of  s»  or  by  any  other  -function  of  s.  The  interpolation  formula  used  here  is 
>jst  a  guesst  or  a  speculation  guided  by  the  physical  understanding  of  coherent 
quantum  1/f  noise  as  a  collective-field  effect,  and  of  conventional  quantum  1/f 
noise  as  an  effect  which  is  not  based  on  the  collective  field  state  of  the 
particles,  but  arises  from  the  individual  field  of  each  carrier. 

The  most  interesting  component  of  the  recombination  current  is  the  surface 
recombination  current  which  plays  a  major  role  in  the  case  of  infrared  detectors 
with  pn  junctions,  in  the  case  of  MIS  devices  this  role  is  not  so  important,  as  our 
calculation  shows.  Nevertheless,  one  should  try  to  reduce  both  the  concentration 
of  recombination  centers  and  the  value  of  the  surface  potential  jump  U.  This  can 
be  accomplished  with  careful  surface  treatment,  and  with  a  good  passivation 
layer.  SiO-  layers  have  been  successfully  used  by  Radford  and  Jones  in 
ion-implanted  and  double  -  layer  epitaxial  HgCdTe  photodiodes13. 

In  general  the  larger  life  time  of  the  carriers  in  M15  devices,  compared  to 
junction  devices  is  due  to  the  absence  of  the  damage  inflicted  by  ion  - 
implantation,  or  by  the  heavy  doping  required  in  double  -  layer  epitaxial 
photodiodes.  The  quantum  i/f  noise  is  invesely  proportional  to  this  life  time. 
Therefore,  MIS  devices  should  have  lower  i/f  noise.  On  the  other  hand,  i/f  noise 
present  in  the  applied  gate  voltage,  in  the  timing  of  the  readout  and  the  value  of 
the  readout  potential  will  be  added  as  a  i/f  noise  source,  if  it  is  present.  In  the 
present  calculation,  however,  this  noise  source  has  not  been  included. 

Any  reoucuon  in  the  concentration  of  tunneling  centers  present  in  the  bandgap 
will  have  a  positive  effect  on  quantum  i/f  noise.  As  we  have  seen  in  Sec.  II,  p  - 
type  devices  should  yield  less  tunneling  via  bandgap  centers.  The  effective  mass 
present  in  the  denominator  of  the  quantum  i/f  noise  formula  in  this  case  should 
just  be  the  effective  mass  of  the  carriers  after  the  tunneling  process,  i.e.,  the 
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effective  mass  of  the  outgoing  carriers  emerging  -from  the  process  we  have 
considered)  or  the  effective  mass  of  the  carriers  coming  in  to  the  process  of 
tunneling  toward  the  centers  in  the  bandgap.  Here  we  have  considered  the 
tunneling  process  as  the  slower  process  which  actually  controls  the  rate  of 
tunneling  via  bandgap  states.  The  capture  of  carriers  by  the  bandgap  states  is 
the  second  part  of  this  compound  process  and  has  been  considered  fast  enough,  sc 
that  it  does  not  limit  the  rate  of  the  total  process.  In  general,  however,  both 
parts  of  the  process  have  to  be  considered  as  a  limitation  on  the  rate,  and  in  this 
case  our  noise  formulae  have  to  be  revised  through  the  inclusion  of  an  additional 
term  similar  to  tne  recombination  noise  term. 

In  the  case  of  very  small  MIS  devices,  where  only  conventional  quantum  i/f 
noise  should  be  present,  we  may  find  lower  noise  in  the  n  -  type  devices,  whose 
bulk  minority  carriers  are  holes  with  much  larger  effective  masses  than  the 
electrons.  This  may  happen  in  spite  of  the  larger  tunneling  via  bandaaD  centers 
located  right  under  the  surface  of  these  devices. 

Finally,  we  would  like  to  empahsize  that  the  present  study  has  attempted  to 
explain  the  basic  concepts  of  quantum  i/f  noise  and  to  illustrate  their  application 
to  MIS  infrared  detectors.  Although  we  have  tried  to  pursue  the  calculation  all 
the  way  to  the  evaluation  of  the  detectivity,  the  data  which  we  used  in  the 
calculation  may  not  be  applicable  in  the  practical  case  at  hand,  and  may  have  to  be 
replaced  with  pertinent  data  in  any  concrete  case. 

The  author  would  like  to  acknowledge  the  help  of  M.  Belasco,  M.  Kinch,  E.  Kelso 
and  R.  Balcerak  in  many  discussions  on  MIS  devices  and  their  noise  problems. 
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VII.  APPENDIX  B 

BACKGROUND  OF  THE  G  -  R  i/F  NOISE  DERIVATION 

At  -first  sight  it  looks  as  if  (1)  is  incorrect  since  (R  -  G)  should  be  proportional  to 
the  trap  density  Nj.  As  will  be  shown  in  the  following ,  the  trap  density  is 
incorporated  into  the  time  constants Tp0  and  Tno. 

According  to  Warner  and  Grung  C i-4,  eq.  <4  -  191)3 

R  -  G  *  vpApvnAnNt  (pn  -  n^/v^ln  +  nt)  +  vpAp(p  +  nt)  (B  1) 

Here,  vp  is  the  average  speed  of  the  holes,  vp  the  average  speed  for  electrons, 
Ap  the  trap  cross  section  for  holes,  and  An  the  trap  cross  section  for  electrons. 
Multiplying  denominator  and  numerator  in  (B  1)  by  N^  and  then  dividing  both  by 

vpApNtvnAnNt  Xields 


R  -  G  *  (pn  -  rr^/nn  +  nj)Tp0  +  (p  +  n^t^] 


(3  2) 


which  15  just  Eq.  (3.23).  in  addition 


T  pQ  =  i/VpApN^  =  1/Cpi  Cp  =  VpApN^^  (B  2a) 

^no  =  i/vnAnNt  =  l/Cn*  Cn  =  vnAnNt;  (B  2b)  . 

as  had  to  be  proved.  We  see  thatTp0  and  ZnQ  are  reai  -  time  constants 
(dimension  :  seconds)  and  that  Cp  and  Cn  are  real  rates  (dimension  :  s-1). 

Note,  however,  that  T=  Tpo  +  ZnD  is  not  the  time  constant  o-f  the  g  -  r 
process.  According  to  Van  Vliet  C15  3  the  true  g  -  r  time  constant  is 


=  C(n  +  nj)vnAn  +  (p  +  n^VpApD  =  N^/C(n  +  n^Cp,  +  (p  +  n^CpD  (B  3) 
It  we  put  l/Tp  =  <n  +  nj)Cn/N^,  l/2"p  =  <p  +  njJCp/N^,  then 


t'  .  ~  ^ 

^  t  ~  +  cn1‘ 


(3  3  a) 


Our  Eq.  (3.23)  is  thus  tully  correct.  It  Nt  goes  to  zero, Tno  and  T po  go  to 
intinity  and  hence  (R  -  G)  goes  to  zero.  Note,  however,  that  T is  independent  or 
N^,  except  tor  a  possible  dependence  ot  n  and  p  upon  N^. 

We  also  note  that  Ap  and  An  are  practically  equal.  For  the  case  that  the  trap 
is  neaatively  charged  when  occupied  by  an  electron,  the  excited  electron  interacts 
with  a  neutral  trap  and  hence  has  a  cross  section  ot  the  order  ot  But 

when  a  hole  is  created  by  having  an  electron  trom  the  valence  band  recombine  with 
the  empty  trap,  the  electron  also  interacts  with  a  neutral  trap  and  hence  has 
practically  the  same  cross  section.  It  the  trap  is  neutral  while  occupied  by  an 
electron,  the  excited  electron  interacts  with  a  positively  charged  trap  and  hence 
has  a  cross  section  ot  the  order  ot  10"^  ~  iO^^cm*-.  But  the  same  is  t~ue  when 


an  electron  from  the  valence  band  interacts  with  a  positively  charged  trap. 
Hence.  Ap  and  An  are  always  nearly  equal. 

However,  vp  is  inversely  proportional  to  imp  and  vn  is  inversely 
proportional  to  so  thatTno  is  proportional  to  and ”2Tpo  to 

(mp*)^^.  Consequently 

^po2/<rpo  +  7:no)2-  mp*/c<mn*)1/2  +  <mp*,1/23 

^no2/<rpo  ^no’2-  mn*/C(mn*)l/2  +  <mp*)1/2Di/2.  (B  4) 

We  need  this  for  the  calculation  of°^  ^  9  ~  p  noise  theory. 


There  are  two  forms  of  quantum  1/f  noise.  In  the  first  place  Cp  and  Cn 
fluctuate  because  Ap  and  An  fluctuate  independently  in  a  i/f  fashion,  according 
to  quantum  I/f  noise  theory.  This  would  yield  Hooge  parameters  d'Hp  and 
where 


Jt/Hp  =  (4^/3;,  )(4Vp2/c2)  =  <4^/37T)(3KT/mp*c2) 

ot'Hn  =  <4^/37I)(/ivn2/c2)  =  (4^-/3?!)(3kT/mn#c2},  (3  4a> 

since  mp*avp2/2  =  mn*zwn2/2  =  3/2KT  for  the  excitation  of  electrons  and  holes 
from  traps. 

In  the  second  place  the  electrons  and  holes  in  the  space  charge  region  are 
decelerated  or  accelerated.  The  potential  across  the  space  charge  region  of  width 
w  eff  *s  a<Vjlif  “  where  a  =  1/2  for  a  symmetric  junction  and  aCz'  0.75  fo  an  n+  - 
p  junction.  Electrons  excited  from  the  trap  level  are  accelerated  and  thus  gain  an 
energy  En  =  ea(Vdi^  -  V)  =  mn*Avn2/2  and  holes  excited  from  the  trap  level  are 
accelerated  and  gain  an  energy  Ep  =  e<i  -  aiUV^  -  V)  =  This 

produces  quantum  1/f  noise  such  that 


(B  5) 


/Hp  =  C2e(l  -  a)(Vdif  -  V)/mp*c23(4ot/3;p 
cl"Hn  =  C2ea(Vdif  -  V)/mn#c2](4^/37T)  (B  5a) 

Adding  the  two  independent  contributions  too£-^n  and^^p  yields 
d  Hn  =  <4«(/37rKC2ea(Vdi<:  -  V  +  3kT3/mn*c2> 
ctHp  =  <4^/37TKC2e<l  -  a)(Vdif  -  V)  +  3kT]/mn*c2},  (B  6) 
so  that  -from  (B  4),  (B  6),  and  (3.35) 

“  ^"po^po  ’’"Z'no^^Hp  +  ^no^^po  +  ^no^-^Hn 
=  (V/37T)(C2e(VdH  -  V)  +  6KT3/C(mn#)1/2  +  (mp#)1/232c2.  (B  6a) 

For  Cdj.Te  with  x  =  0.20,  mn*  =  0.01m  and  mp*  =  0.55m  so  that  C(mn*)1/2  + 
(mp#)1/2]2  *  0.7im  which  is  a  -factor  7i  larger  than  m*n.  An  earlier  estimate  by 
Radford  and  Jones  C 1 3 D  used  mn*  instead  of  E(mn*)1/'2  +  (mp*)^2]2. 
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